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Abstract 



This thesis is concerned with problems related to Synthetic Aperture Radar 
(SAR), a technique of making images of the surfaces of planets using electro- 
magnetic waves. Reconstructing images of the surfaces from the gathered data is 
an inverse problem as are other young and thriving imaging techniques used for 
example in optical tomography and transient elastography. In optical tomography 
one tries to create images of the human body employing light, whereas in tran- 
sient elastography ultrasound is used to measure the propagation of shear waves 
and thus reconstruct the stiffness of human tissue. But the field of inverse prob- 
lems covers older and established topics as well like computerized tomography 
that creates images of the human body by means of x-rays and magnetic reso- 
nance imaging using electromagnetic fields to measure the distribution of atoms. 
All these very different applications have in common that the gathered data is 
difficult to interpret. Therefore mathematical processing is necessary in order to 
create an intelligible image of the measured object. Some of the problems related 
to this mathematical processing necessary in creating SAR-images are analyzed 
in this thesis. 

The thesis is structured as follows: The first chapter explains what SAR is, and 
the physical and mathematical background is illuminated. 

The following chapter points out a problem with a divergent integral in a common 
approach and proposes an improvement. Some numerical comparisons are shown 
that indicate that the improvements allow for a superior image quality. 
Thereafter an important problem is analyzed - the problem of limited data. In 
a realistic SAR-measurement the data gathered from the electromagnetic waves 
reflected from the surface can only be collected from a limited area. However the 
reconstruction formula requires data from an infinite distance. The chapter gives 
a comprehensive analysis of the artifacts which can obscure the reconstructed 
images due to this problem. Additionally, some numerical examples are shown 
that point to the severity of the problem. 

In chapter |4] the fact that data is available only from a limited area is used to 
propose a new inversion formula. This inversion formula has the potential to make 
it easier to suppress artifacts due to limited data and, depending on the application, 
can be refined to a fast reconstruction formula. 

In the penultimate chapter a solution to the problem of left-right ambiguity is 
presented. This problem exists since the invention of SAR and is caused by the 
geometry of the measurements. This leads to the fact that only symmetric im- 



ages can be obtained. With the solution from this chapter it is possible to recon- 
struct not only the even part of the reflectivity function, but also the odd part, thus 
making it possible to reconstruct asymmetric images. Numerical simulations are 
shown to demonstrate that this solution is not affected by stability problems as 
other approaches have been. 

The final chapter lists some conclusions drawn from the preceding chapters and 
develops some continuative ideas that could be pursued in the future. 

Acknowledgements 

First, I would like to thank Prof. Dr. Dr. h. c. F. Natterer for stimulating the work 
on this thesis and for offering helpful advice. 

My thank also goes to Dr. F. Wiibbeling who contributed to this thesis in many 
fruitful discussions and to all members of the institute for the pleasant atmosphere. 
Finally I would like to especially thank my parents and my significant other, 
Birgit, for their ceaseless support that made this thesis possible. 



Contents 



1 Introduction 

1 . 1 Inverse problems 

1 .2 Physical background 

1 .3 Mathematical model 

2 An improved inversion formula for the spherical Radon transform 

2.1 Introduction 

2.2 Properties of the function g = Rf 

2.3 Modification of Andersson's first inversion formula 

2.3.1 Definition and properties of a modified R* 

2.3.2 A modified inversion formula 

2.4 Problems with Andersson's second inversion formula 

2.5 Numerical simulations 

2.5.1 New approach 

2.5.2 Comparison 

3 Ghosts due to limited data 

3 . 1 The problem of limited data 

3.2 Orthogonal functions and their transforms 

3.3 Projection of unmeasurable data onto orthogonal functions . . . 

3.4 Reconstruction of the orthogonal functions 

3.5 Numerical simulations 

3.5.1 Reconstructions for g e {o^f : a 6 R, b > 0} 

3.5.2 Reconstructions for g e {o^ en : a > 0, / G N } 

4 A new approach to invert the spherical Radon transform 

4. 1 Orthogonal functions and their transforms 

4.2 Projection of the data onto the orthogonal functions 

4.3 Reconstruction of the orthogonal functions 



I 



5 Tackling the left-right-ambiguity 

5.1 Using two data sets 

5.2 Post-processing formulas 

5.3 Numerical simulations 

5.3.1 Solving the left-right ambiguity with computed symmet- 
ric images 

5.3.2 Solving the left-right ambiguity with reconstructed sym- 
metric images 

5.3.3 Analysis of stability and the effects of noise 

6 Conclusions and outlook 
Bibliography 





73 




74 




75 




78 




79 




82 




87 




99 


101 



II 



Chapter 1 
Introduction 



This chapter will give a short introduction into the mathematical and physical 
concepts necessary to deal with SAR. After a closer look at the mathematical 
terms of well-posed versus ill-posed and direct versus inverse problem, the second 
section will describe the physical and technical background related to SAR. The 
last section will then give an introduction to the mathematical model of the inverse 
problem involved in SAR. 

1.1 Inverse problems 

From a mathematician's point of view the problems associated with SAR be- 
long to the field of inverse problems. It is difficult to exactly define the term 
inverse problem. Therefore to begin with, the closely related ideas of well- and 
ill-posedness and the concept of stability will be illuminated. To this end a trip 
into history might prove helpful. The scientific community was aware of the prob- 
lem of instability as early as 1873 when Maxwell wrote in an essay from February 
11th, 1873 dp. 434]: ' 

"There are certain classes of phenomena, as I have said, in which a small error 
in the data only introduces a small error in the result. ... The course of events 
in these cases is stable. There are other classes of phenomena which are more 
complicated, and in which cases of instability may occur, the number of such 
cases increasing, in an exceedingly rapid manner, as the number of variables in- 
creases. " 

The concept of a well-posed problem was first formulated by Hadamard in 1902 
[0. A well-posed problem according to Hadamard requires the existence, unique- 
ness and stability of a solution, originating from the philosophy that the mathe- 
matical model of a physical problem has to have these properties. If one of the 
properties fails to hold, a problem is called ill-posed. Due to this background only 
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well-posed problems were studied extensively for quite some time and ill-posed 
problems were neglected. 

This led to a vast knowledge of and familiarity with direct problems which are 
usually well-posed. These are problems that are common in physics with an input 
x and an operator H modelling the physical system. In these cases it is the goal 
to find the system's response y = Hx. Therefore they are called direct problems 
as it is the classical physical problem to have a system state and let it evolve ac- 
cording to a certain formula. The term inverse problem is derived from this notion 
of direct problems. The corresponding inverse problems consist in determining 
the cause x for a known system H and a known response y or in determining pa- 
rameters of the system H for a known input x and a known response y. Unlike 
direct problems that are usually well-posed, inverse problems are often ill-posed. 
Inverse problems are relevant for many different physical applications, for exam- 
ple problems related to imaging. One of these applications is SAR which will be 
described in detail in the following. 

1.2 Physical background 

Synthetic Aperture Radar is a technique of taking pictures of the surface of planets 
from an airplane or satellite. As the name indicates, SAR utilizes Radar, electro- 
magnetic waves with a much longer wavelength than used in optical imaging. 
These waves are emitted by antennas mounted on the airplane or satellite. They 
are then reflected from the surface and detected by the same antennas. The use 
of such large wavelength leads to a great advantage in comparison with photog- 
raphy. Since waves of these wavelengths can penetrate clouds and even foliage, 
SAR images can be taken in foggy or cloudy weather. This capability is demon- 
strated in figure [TTI that shows in both images the region of Waterford in Ireland 
on the morning of August 9th, 1991. The left picture shows an optical image that 
is almost completely obscured by clouds whereas the SAR image on the right is 
not at all affected by them. This is an important feature since in Europe only one 
out of ten optical images is free of clouds [3]. However, this comes at the cost of 
a reduced resolution compared with optical imaging. 

The typical frequency used in SAR is largely dependent on the technical imple- 
mentation. Commonly frequencies ranging from 20 MHz to 10 GHz are used [4], 
0. This corresponds to wavelengths of 15 m to 3 cm. They even permit to detect 
concealed object, e. g. covered by trees, and to measure the biomass of a region 
0. 

An example of a SAR system is CARABAS. This is an airborne VHF SAR system 
developed by the FOA (National Defence Research Establishment) in Sweden. 
CARABAS consists of two 5.5 m long antennas mounted parallel on a Rockwell 
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Figure 1.1: Waterford (Ireland) on the morning of August 9th, 1991. Left: Optical 
Landsat satellite image. Right: SARERS-1 satellite image. Source: 



Sabreliner aircraft. The distance of the antennas corresponds to the shortest emit- 
ted wavelength and translates to three meters since the used frequency spectrum 
ranges from 20 to 90 MHz. The airplane usually travels at a height of 1500 to 
10000 m at a speed of 100 to 130™ 01. For a frequency of 70 MHz theoretically 
the best achievable resolution is 1 m parallel to the flight track and 2 m perpendic- 
ular to the flight track [6]|. However it turned out that for real measurements the 
resolution obtainable is only half as good Q. 

1.3 Mathematical model 

The correct mathematical model for radar emission and scattering is given by 
Maxwell's equations of electromagnetism. Therefore Maxwell's equations are 
used to understand certain effects that are common in SAR and to enhance the 
reconstruction IfTOl . But due to the complexity of deriving a reconstruction 
formula based directly on these vector equations, usually only a scalar wave 
model is used in mathematical models (H, flU. In today's real life applications 
not even algorithms based on the wave equation are commonly used, but a simple 
backprojection algorithm is applied that is now explained in more detail. 

For the simplest case assume a point source that is flying along a straight line 
above a flat plain and emitting infinitely short pulses in constant time intervals. 
For a single pulse this results in an expanding sphere centered at the emitting 
position. When the expanding sphere is large enough, it hits the ground and is 
reflected. Since the speed of light in air is much larger than the speed of the 
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Figure 1.2: The geometry of SAR. Source: IfTTTl 



airplane, the airplane is considered static during each send-receive cycle (start- 
stop approximation). As the scattered wave is therefore received at exactly the 
same position at which it was emitted, the received signal comprises an integral 
of the ground reflectivity over the circle that is obtained intersecting the plane 
and the expanding sphere at a certain point in time as depicted in figure 1 1 .21 The 
equation that describes this measurement is 



where g is the measured data and / is the function to be reconstructed (this no- 
tation will also be adopted in the following chapters). S 1 denotes the unit sphere 
in R 2 , dSi stands for the canonical surface measure on 5 1 , r denotes the radius 
of the circle of the intersection of the plane and the expanding sphere, and x is 
the position of the airplane on the track projected onto the plane. This equation is 
called the spherical Radon transform. It is obvious that the collected data contains 
only the even part of the reflectivity function with respect to the flight track as the 
whole setup is left-right symmetric. The inverse problem now consists of recover- 
ing /, either analytically by finding a reconstruction formula or numerically. The 
easiest way to do so is a simple backprojection which numerically distributes each 
point of collected data equally over all points of the reconstructed image that lie 
on the circle where the data comes from. 

The first thorough mathematical treatment based on this model was developed by 
Andersson lfT2l . Chapter[2]is based solely on these ideas and analyzes them more 
thoroughly. This model is also the foundation of chapters [3] and HI but chapter [3] 
should also yield insight into a more general problem associated with SAR and the 
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ideas of chapter 0] may be transferred onto other approaches to SAR. In chapter[5] 
the algorithms based on this simple model are only used for numerical examples, 
whereas the underlying theory developed in this chapter is much more general. 
Therefore the ideas of this chapter should be applicable to all SAR systems. 
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Chapter 2 



An improved inversion formula for 
the spherical Radon transform 

The determination of a function from spherical averages is a problem often en- 
countered in physical applications such as SAR and SONAR (SOund Navigation 
And Ranging is a technique that uses sound propagation under water to navigate 
or to detect other watercraft). The work related to this topic, which has lead 
to a great amount of insight and refinement today, began with the proposal of 
a reconstruction formula by Fawcett lfP3l . The mathematical analysis of the 
problem was later improved by Andersson Ifl2ll . and two refined reconstruction 
formulas were derived. This sparked a host of activity |[T4ll . 0, lfT5ll . IfToll . IfTTll 
so that much of today's research is based on Anderssons's ideas. 

In lfT2ll two reconstruction formulas were derived from the Fourier inversion for- 
mula, but it was not checked whether they are properly defined. In the following it 
is shown that the first contains an integral that diverges under physically sensible 
conditions. An alternative is presented. Additionally it is shown that the other 
reconstruction formula might be difficult to compute numerically. 

2.1 Introduction 

At first, for the benefit of the reader, some results from lfT2l will be recalled. Note 
that besides the aforementioned problems, a few minor errors occurred, which do 
not essentially obscure the results in Ifl2l . For a detailed analysis of these errors 
see lfT8l . Additionally it will be shown that an integral in the first reconstruction 
formula in Ifl2ll does not converge. 

2.1.1 Definition 

LetneN. Then 
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1. N := NU {0} 

2. ,y(R n ) := {if G C°°(R n ) : ||x /3 d t VlU < oo V/3, a G N™} is the Schwartz 
Space. 

3. A sequence (<^z)zeN C ^(R n ) is caiied zero sequence if 

Wx^ipiWoo -> 0V/3,a G Nq. 

Z— >-oo 

4. A linear functional f : ^(R™) — >■ C, </? h-K /,</?> is continuous if 
< f,(fi > — ► for all zero sequences (<^z)z s n C ^(R n ). 

Z— »oo 

5. J^ e (R n+1 ) := G J^(R n+1 ) : <p(x,-y) = <^(x,y),Vx E R™, y E R}. 

6. ^ r (R n x R n+1 ) := G ^(R 2n+1 ) : V orthonormal transformations 

U : R n+1 -)• R rt+1 Vx eR",2 6 R n+1 : ^(x, z) = ?7z)}. 

7. J^'(R n+1 ) and J^ r '(R n x R n+1 ) are the duai spaces of ^ e (R n+1 ) and 
^ r (R n x R n+1 ), respectively. The weak-* topology is used. 

8. (y;)„ e (R"xR" +1 ) := {g E ^(R"xR" +1 ) : supp g C {(£,n) : ||n|| > 
lleil}}- 

9. S 1 ™ denotes the unit sphere in R" +1 , and S n denotes the canonical measure 
on the surface of S n . 

10. Let f E ^(R n ). Then f denotes the Fourier transform in R n : 

R" 

11. Let f E y(R n ) and g E <¥"{W). Then the Fourier transform of g is 
defined by 

<9,f >=<9,f > ■ 

12. Let f E y e {R n+1 ). Then the operator R is defined by 

Rf(x, r) := J f(x + r£, rr]) dS n (£, 77). 
8 



In this definition the operator R describes, in the case of SAR and SONAR, the 
measurement of the reflectivity function / that represents the ground reflectivity. 
The measurement is modeled as a 5-impulse wave that propagates as concentric 
spheres. The ground is approximated as a plane. The single scatter approximation 
for the wave hitting the ground results in integrals over circles. 
For simplicity, in the following sometimes only y, S" e , etc. is written instead of 
y(R n+1 ), y e (R n+1 ), etc. 

2.1.2 Remark 

1. Let f G y e (R n+1 ). Then for g = Rf G J^'(R n x R n+1 ), x G R n , and 
r > in the following g(x, r) is sometimes written with abuse of notation. 
This identification ofg(x, y) with g(x, \y\) fory G R n+1 is justified because 
g depends only radially on the last n + 1 variables. 

2. Since ,y r (R n x R n+1 ) and y e (R n+1 ) are subsets of J^(R n x R n+1 ) and 
y(R n+1 ), the definition of the Fourier transform can be extended to y r , 

3. It is easily seen that f G y e andg G y r iff G y e andg G y r , respectively. 
The essential result in lfT2ll is the (Fourier-) inversion formula: 

2.1.3 Theorem 

If y e (R n+1 ) is given the topology of y'(R n+1 ), the mapping 

r : y e (R n+1 ) y;{u n x n n+1 ) 

is continuous and can, by continuity, be extended to a mapping 

R : y^(R n+1 ) -> y' r {R n x R n+1 ). 
The range of this extended mapping R is the closed subspace 

(y;) cone (R n x R n+1 ) c y r {R n x R n+1 ). 

R is one-to-one and the inverse mapping 

r- 1 ■■ {yXone{R n x R n+1 ) y^(R n+1 ) 

is continuous. Moreover, if g = Rf and /(£, rf) or g(£, rf) are integrable functions 
for £,7/ in some open set, then in that set 



2 / e,VNl 2 -lieiP y 

forO<\\r)\\<U\\ 



or 

1 \Qn\ „_! . x 

= y m (lien 2 + ^ $ (e, vIm) . 

Proof: 

lfT2l Theorem 2.1] 

□ 

This means that the Fourier transform g of the data g can under certain conditions 
be used to extract the Fourier-transform / of the reflectivity function /. 

2.1.4 Definition 

1. Letg e y r {R n x R n+1 ). Then 




R" 



2. Let f G L 2 (R). Then the Hilbert transform of f is defined by the principal 
value integral 

ir J x-y 

R 

In lfl2ll two reconstruction formulas were derived from theorem 12.1 .31 the first 
of which was already essentially given in [fT3l . They are given in the following 
corollary. 

2.1.5 Corollary 

1 \S n I 

With c n = -^y^ l y ± two reformulations of the inversion formula are possible: 

1. Forg G y r 

f = c n H y -^-A^R*g 

with the Hilbert transform in y, H y , and the Laplace-Operator A = A x + 
^2 . This formula is essentially also given by Fawcett HJ\l . 

2. Forg e (^) co „. e := {9 <= ,9> r : supp g C {(£,77) : > U\\}} 

f = c n R*Kg 

with the operator K defined by Kg(£, rj) = \/IMI 2 ~~ ll£ll 2 IMI™'~ 1 #(£; V)- 



10 



Proof: 

Ifl2l Section 3] 

□ 

The corollary states that under certain restrictions it is possible to reconstruct the 
reflectivity function / directly from the data g without taking the detour through 
the Fourier space. 

2.1.6 Remark 

Note, that the essential restriction of this formulation of the corollary in compar- 
ison to [T2\l is that the data g has to be in y r or (y r )cone respectively. This is 
necessary because otherwise the application of R* to g in the first case is unde- 
fined or the application of K to g in the second case. 

Unfortunately g = Rf is usually not in 5^ r . Therefore the two reformulations of 
the inversion formula are only valid in the distributional sense with an appropri- 
ately defined R*. 

For example a physically reasonable f £ y e , f : R 2 — > [0, oo) with f(x ,yo) > 
for some x , y £ R yields f(x, y) > c > for all (x, y) £ K € (x , y ) with ap- 
propriate c, e > 0. 
Moreover 

2ng (z, y/(z - x) 2 + y 2 ) = [ f( y z + r 1 ,r 2 )da(r) 
v ' \l{z — x) lJ ry l J 

\\r\\ = ^{z-xYW 

with r = (ri, r 2 ) and a the canonical measure on the sphere with radius r in R 2 . 
For \/ (z — x ) 2 + > e we obtain with a simple geometrical consideration: 



2ixg (z, \J(z- x ) 2 + 



ec 

> 



Thus 



27rR*g(x ,y ) = / g I z, ^ (z - x ) 2 + yl J dz 



> / — dz = oo. 



2 



^/{z-x ) 2 +yl >e 



\/{z-x ) 2 + yl 



Since this integral diverges, the reconstruction of a non-negative function f with 
f(xo, yo) > for some (x , yo) £ R 2 is impossible using the first reconstruction 
formula. 
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This result is in accordance with a result from Nessibi, Rachdi, and Trimeche lfT9l . 
They gave reconstruction formulas for functions g = Rf with 

J P(y)f(x,y) dy = 
o 

for all x G R n and for all one-variable polynomials P. 



2.2 Properties of the function g = Rf 

Before showing an important property of the data function g = Rf, which will 
be necessary for the derivation of the new reconstruction formula representing an 
alternative to Andersson's, some definitions will be needed. 

2.2.1 Definition 

Letn G N. Then 

1. 3l(R n+l ) := C^(R n+1 ). 

2. 2$ e (R n+1 ) := {cp G ®{R n+1 ) : <p(x,-y) = ip(x,y),Vx G R™, y G R}. 

2.2.2 Lemma 

Iff G y e {R n+1 ), then g = Rf G C°°. 
Proof: 

The interchange of differentiation and integration is justified because / G 

y e (R n+1 ). 

□ 

Now an important property of g can be shown. 

2.2.3 Theorem 

Iff G y e (R n+1 ), then g = Rf G L 1 (R n x R n+1 ). 
Proof: 

With fey, also / ey C L 1 . Theorem [2"T3l implies that 

2 /(wNi 2 -iieii 2 ) 

ut n ) - ) (2vr) n -— -4 7 for V > £ 

\ Sn \ Nh~VIMI 2 -ll£ll 2 

otherwise. 



Therefore 



|»(e,»7)|de^ = (27r) B 1 2 



IS* I 

R"xE»+ 1 INI>II£II 
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/ U,V\\v\\ 2 -UW- 



dl; drj. 



The substitution \\r)\\ = p' results in 



The substitution p' = p + ||£|| leads to 



p'fU, Vp' 2 -UW- 



Vp' 2 -UW' 



d£ dp'. 



m,r ] )\dZdr ] = 2(2nr 




R"xE" +1 

/ is in 5? , therefore 

\g(£,ri)\d£dri 



< 2(2tt) 



R" p>0 



(p+U\\)fU,y//^ + 2 P U\\ 



v / p 2 + 2p|iei 



d^ dp. 



R"xE' 1 + 1 




c(p+U\\) 



R" p>0 



2(2tt) ? 




v/p^ + 2 p|K|| (l + v/|K|| 2 + p 2 + 2p|K|| 
C(p+||ell) 



n+3 



d£dp 



R" p>0 



<2(2nrJ J 



v V + 2p||£||(l + ||eil+p)" +3 

c 



R" p>0 



v^+Wlfa + lieii+p) 



n+2 



and the integral 




C 



R n p>l 

converges. So only the integral 



v /p 2 + 2p|ieii(i + neii + p)" 



+2 



d£ dp 



C 



R™ 0<p<l 

remains to be examined. 



R" 0<p<l 
< 



v V + 2p||£||(l + ||£||+p)"+ 2 

C 

vV + 2p||£||(l + ||£||+p)«+ 2 

C 



d£ dp 



d£ dp 



R n 0<p<l 



'pum+uw) 



n+2 



d£ dp 
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/ 



v^(i + lieil) n+2 



d£ < oo. 



□ 



2.2.4 Remark 

An analogous proof shows r\g £ L 1 . 

2.3 Modification of Andersson's first inversion for- 
mula 

To derive an inversion formula that overcomes the problem of the diverging inte- 
gral, a modified operator R* is defined. This enables a convenient formulation of 
a new inversion formula. 

2.3.1 Definition and properties of a modified R* 

Now a modified version of the operator R* is introduced, and its properties are 
discussed. 

2.3.1 Definition 

For / £ @ e and g = Rf £ C°°(R n x R n+1 ) we define 



This slight modification by an additional derivation turns out to ensure the conver- 
gence of the integral applied by the operator R* d under the minor and physically 
feasible constraint that / is in 3) e . Thereby the formulation of a mathematically 
exact reconstruction formula is possible. 

2.3.2 Proposition 

R* a g is well defined for f £ St e and g = Rf. 
Proof: 

Let x £ R n , y £ R. Then 





x - z\\ 2 + y 2 ,r]>/\\ 



x - z\\ 2 + y 2 ^j dS n (£,r)) dz . 
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/ G 3) e and therefore it is possible to interchange differentiation and integration. 

\(R* d g)(x,y)\ 



S»\J v/||a;-z|| 2 + r 

R" S n 



< 



x r ) ■ (v/) (z + Zy/\\ x - z \\2 + y2, n^Wx-zr + y*) dS n (£, v ) dz 



R" S™ V 



X 



(V/) (z + Zy/\\x-z\\* + y*,riy/\\x-z\\z + y^j dS n (Z,ri)dz. 



The substitution r = (£, 77) a/J[x — z\\ 2 + y 2 yields 
\(R*a9)(x,y)\ 



< 

TR 



1 



5* 



||r|| = v /||x-2|P+ y 2 

Since f E 

\(R*a9)(x,y)\ 



^/\\ x -z\\ 2 + y^ 



^||(V/)((z,0) + r)||da(r)d*. 



< / z: x ^+1 max (II Wll) diam(supp /) cfe < oo. 

□ 

2.3.3 Corollary 

For f £ @ e andg = Rf, R* 9 g G C°° H 

Proof: 

This is guaranteed by the estimate in proposition l2.3.2l 

□ 

2.3.4 Definition 

Letf G y(R n ) and if e J^(R n ). Then 

< >y(u")-= f(<p) 

is the functional f applied to the test function if. Here the subscript y(H n ) is a 
reminder of the space of the test function ip. 
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Now in analogy to lfT2ll an expression for R* a g is derived. 

2.3.5 Theorem 

Let fe%andg = Rf. Then R* a g(£, rj) = irjg (f , vfpT? 

Proof: 

Let ip E y, f E @ e , and g = Rf. Then 

< R* 9 g,¥ >y(rxE)=< R*a9,V >,/(K»xR)= J {R* d g){x',y)0(x',y) dx' dy 



R"xR 



= / J ^9 (-.VW^WTf) dz0(x',y)dx'dy. 

E"xRR" 

As dp E y and with corollary 12.3 .31 Fubini's theorem implies 

<W>, (M) =/ / (c, VW^W^f) 0(x', y) dx' dy dz. 



R"R"xR 

Continuing in the distributional sense 







< R* d 9, V >y(R»xR)= - / / 9 (z, \/\\x' - z\\ 2 + y 2 ^ -k-<p(x', y) dx' dy dz 



dy 



R™ R" X R 




g (z, V\\ x ' - z \\ 2 + y 2 



R™ R"xR 



d_ 

dy I (2tt 



f ^ <x '^ + ^^)didrAdx'dydz 

) 2 J \ 

R"xR / 




(2tt) 

R" R n X R R" X R 



_ _ e -i<z£> e -i(<x'-z,Z>+y V ) 



x rig [z, \/\\x' — z\\ 2 + y 2 ^ r]) d£ drj dx' dy dz. 
The substitution x' = x + z yields 

< Rfg ,v> >y(rxR) = i J J J ^- ¥& -^>^<^>+w) 



(2tt) 

R" R"xRR"xR 



x rig (z, a/||x|| 2 + y 2 j (p(£, 77) d£ dr\ dx dy dz 
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R-xR 

=< irjguf >/(e»xr) • 

Here g\{£, rj) := g U, y/\\£\\ 2 + V 2 ) and with this definition g\ G ,9" {R n x R). 

□ 

2.3.2 A modified inversion formula 

With the results of the preceding sections a well defined reconstruction formula is 
also attainable for physically meaningful reflectivity functions. 

2.3.6 Theorem 

n — 1 

Let f G and g = Rf. Then f = c n H y A~2~ R* d g with the constant 
c n := j2^~2 an d the Hilbert transform H y . 

Proof: 

Let / G @ e . It follows from theorem [2X3] and theorem [2X3] that 

/(£, V) = c n \v\ (lief + rf) V g (e, VU\\ 2 + V 2 ) 

= c n (-i) sgn{rj) (||£|| 2 + r, 2 ) ^ irjg (£, vW + ?) • 

Theorem 12.3.51 yields 

Tl-l _ „ 

m,v)=c n (-i)sgn(r 1 )(U\\ 2 + r J 2 ) 2 R* d g^, V ) 
= c n (-i) sgn(rj) (A^fl^) rj) 

= c n (H y ^R* a gy^r l ). 
This completes the proof, because / is in S^ e . 

□ 

2.4 Problems with Andersson's second inversion 
formula 

In the following it is shown that the second reformulation of the inversion formula 
in lfT2l is only valid in the distributional sense. 
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2.4.1 Proposition 

Let f ey e mdg = Rf. Then Kg e x R n+1 ) for all p>0. 



Proof: 

Theorem 12.1.31 implies that 

f 2 fte.V\M\ 2 -K\\\ 

w M rMlFP for|l " ll> 11511 

I otherwise. 



Hence 



(2vr) n r^/ (e, vfFIF) for IMI > ||£|| 
otherwise 

2ra+2 



<< W^( 1 + vfMNFffl)"" 7 for ||r?|| > Hell 
otherwise 

otherwise 

<(27rr^(l + max(||e||,h||)r^ 

with an appropriate C > 0. 

□ 

This proposition is sufficient for the validity of the second reconstruction formula 
in the distributional sense.However it is noteworthy that in general g is only in C°° 
and not for example in L 2 . Therefore g has to be computed in the distributional 
sense and no further improvement for this inversion formula is achievable. 

2.4.2 Corollary 

Let f e% with /(0, 0) ^ and g = Rf. Then g,g £ L 2 (R n x R n+1 ). 
Proof: 

Theorem 12.1.31 yields 

2 f(tV\M\ 2 -U\\ 2 \ 

«{,,) = \ ' 2 '»M H^inf for """ >m 

) otherwise. 
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Therefore 



m,v)\ 2 dtd V = (2n) 



2n 



The substitution 



\S n \ 



f (e, v\\v\\ 2 - uw- 



2n-2 



\H\>\\Z\\ 



d£ drj. 



p' results in 



\g(t,r))\ 2 dtdr) 



(27T) 



2n 




J2-n 



fu,Vp l2 -uw- 



K"p'>||€|| 

The substitution p' = p + ||£|| leads to 



p' 2 - neiP 



R n x]R n+l 



(27T) 



2d 




(P+IKII) 2 - / U'V / P 2 + 2p||ell 



P 2 + 2p|KII 



<i£ dp. 



We consider 



(p+neii) 2 - n /(e, vV + 2pikii 



P 2 + 2p|KII 



<i£ dp 



and assume without loss of generality 

l, 

2' 



0<p<i: 



/(e, V / P 2 + 2p||e||)| > lfor neii < i 



(p+neii) 2 - n /(e,Vp 2 +2p|ieii 



ll?ll<4o<p<i 



> 



l<|o< P <i 



p 2 + 2p|ieii 

(p + !ieii) 2 ' ra 
p 2 + 2p|ieii 
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d£ dp 



d£ dp. 



With the substitution r = ||£|| we obtain 



(p+U\\) 2 ~ n ffay/f + MZW 



ll«ll<^0<p<| 



> \s 



n-l| 



0<r<i 0<p<± 



p 2 + 2p|ieii 

(p + r) 2 ~~ n r n ~ 1 
p 2 + 2pr 



d^ dp 



dr dp. 



For n = 1 this is 



(p + MW) fU>y/? + M\t\\ 



iieii<|o< P <i 
> 

0<r<| 0<p<§ 

For n > 2 this results in 



P 2 + 2p||e|| 
p + r 



d^ dp 



p 2 + 2pr 



dr dp > 



0<r<l 0<p<§ 



— dr dp = oo. 
2p P 



(p+u\\) 2 - n fu,vp 2 + 2 p\m 



\<ho< P <h 



P 2 + 2 P U\\ 



d£ dp 



> 



„n-l 



0<r<h 0<p<§ 



p 2 + 2pr 



dr dp = oo. 



□ 



2.4.3 Remark 

Unfortunately usually /(0, 0) ^ for a reflectivity function f ^ with physically 
realistic properties. Moreover g is not continuous, and therefore g L 1 , so it 
might be difficult to compute g numerically with sufficient accuracy. Therefore 
the other inversion formula seems to be a better approach numerically. 



2.5 Numerical simulations 

As can be seen in theorem [2.3 .61 the exact reconstruction of the reflectivity func- 
tion / requires data from the whole half plane. This is impossible in practice, 
therefore several ways to handle this problem were developed, ranging from a 
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simple cutoff to applying an exponential decay towards the edges of the data. In 
the following a new approach is proposed to lessen the artifacts caused by this 
limitation of the data. Then the results of computer simulations using this new ap- 
proach are compared to the results using a simple cutoff. The problem of limited 
data in the spherical Radon transform is thoroughly examined in chapter |3j 



2.5.1 New approach 

The proof of proposition 12.3.21 yields 



i 



y 



\S»\J ^|| x _ 2 ||2 +y 2 



(V/) (z + ^\\x-z\\* + y\ V y/\\x-z\\* + y^ dS n (£, 77) dz 



^ / i 

\\ r \\=s/\\ x - z \\ 2 +y 2 



n+1 



■ (V/)((z,0)+r) da(r) dz. 



The new approach presented here is based upon the idea that the integral 



I-.(V/)((z,0) + r) da(r) 

\\r\\=y/\\x-z\\*+V 2 

does not vary much for large values of z because for large z the integral describes 
a circle with a large radius that runs through the support of /. Therefore the 
curvature only changes slightly and since / e the same should hold for the 
integral. The missing data is replaced by an approximation that uses the first and 
last known data with regard to the variable of integration z as an approximation 
for the interval of integration where the data is unknown. Since it is easily seen 
from the formula above that ^g{x,y) = for y = 0, it is sufficient that the 
following approximation is well defined for y ^ 0, where z min and z max denote 
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the smallest and the largest values of z with data available: 

Zmax 

—g (z, y/\\x - z\\ 2 + y 2 ^j dz 



^ m i n 



y 



\s n \ 



y/\\x - z\\ 2 + y 2 ^ 



n+1 



X Zmin \ \ D ^ 



n+1 



y ( . /iu _ z . 112 _l „,2 N n 



X 



I X Zmin || V J 

J \\ r \\ 



\\r\\=-\/\\x—Zmin II 2 +V 2 



y/\\x- z\\ 2 + y 2 ^j 



Zmax 

n+1 



n+1 



\/p - z max \\ 2 + y 2 ^ 



n+1 



\x Zrnax \ | 2 J 

x y (Vf)((z max ,0) + r)da(r)dz 

\\ r \\ = >/\\ x — z max\\ 2 +y 2 

Zmax 



>/P - 2 min || 2 + ?/ 2 j 



n+1 



a/||x- z|| 2 + y 2 ^j 



n+1 



X (zmim \A\ X ~ z min\\ 2 + ?/ 2 ) cfo 



n+1 



(v/||x-^| 2 + y 2 ) 



n+1 



x £-g (z max , vf^JPT?) d* 
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Zmax 




2.5.2 Comparison 

Now the images obtained via the approximate reconstruction formula derived 
above are compared to the results of a simple reconstruction formula that sets 
unmeasured data to 0. This comparison will highlight the advantages of this new 
approach. 



f(0,y) 




Figure 2.1: Phantom 

The underlying phantom for the following comparison is very simple. It is a circle 
with a radius of 20 as depicted on the left in figure I2T1 Its center is with (0, 25) 
close to the flighttrack that runs along the left edge (y = 0). The circle has 
an amplitude of 10, the remaining part has a reflectivity of 0, as can be seen in 
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the cross section on the right that is indicated in the left image by the horizontal 
line. The figures shown in the following are all composed in this same way. The 
simulated measurements are noise free. Several reconstruction pairs are shown. 
The first uses data that stretches exactly as far as the reconstructed area, i. e. 
< r < 256 and —128 < x < 128. The following pairs are reconstructed each 
with a larger amount of data in both directions with respect to the preceding one. 
The reconstructions are designated accordingly. The first reconstruction in each 
pair is obtained by continuing the data with in the region where data is missing. 
The second is reconstructed using the approximation derived above. 
Figures 12.21 and 12.31 show reconstructions that use an amount of data that is ex- 
actly as large as the data contained in the images, i. e. < r < 256 and 
— 128 < x < 128. Data that is unavailable is set to for reconstruction pur- 
poses in figure 12721 This causes two broad circular artifacts - one curved upward 
and one downward - that are clearly visible. The formation of the artifacts can 
be understood as follows. The information in the data that is actually measured 
causes the algorithm to reconstruct the large positive circular phantom. The only 
possibility however to conform to the of the continuation of the data is to form 
these circular negative artifacts. This is the only way to achieve that the circular 
integral which runs through the large positive phantom becomes 0. Therefore the 
radius of these artifacts matches exactly the smallest radius that is missing in the 
data. A side effect is that both algorithms overshoot the amplitude of the circle. 
It is notable that in figure [2721 the error is large close to the object and decreases 
as the artifact closes to the edges of the image. Another problem is the noticeable 
gradient in the object's amplitude. 

The missing data was dealt with as delineated in subsection [23T] to obtain figure 
12.31 Again this causes artifacts that are similar to the artifacts in figure 12.21 The 
reconstruction is very good close to the object, as can be seen in the cross section 
before the dip, but the artifacts get worse as they approach the edges of the im- 
age. This comes from the fact that for large values of |x| the approximation gets 
less reliable. However, the gradient in the object's amplitude in figure [2721 is not 
reflected in figure [2731 where the amplitude is constant, as it should be. 
In figures 12741 and [2751 double the amount of data in each direction is used relative 
to the preceding two images, i. e. < r < 512 and —256 < x < 256. Again the 
first figure shows the reconstruction using a continuation by for missing data, 
whereas the second figure is computed with the help of the approximation. The 
artifacts seen in figures [2741 and 12 .51 are similar to the artifacts seen in the previous 
image pair, albeit less severe. The gradient of the circle in figure 12741 is less steep 
than in figure [2721 and the dip below is shallower but broader. A similar effect is 
noticeable in comparison of figures [2731 and [2751 The area close to the object with 
its very good accuracy between the object and the dip, where the reflectivity stays 
constant around 0, is larger and the dip is more shallow. 
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Figures l231 and l2~77l are computed with quadruple the amount of data relative to the 
former pair, i. e. < r < 1024 and —512 < x < 512. The trend of the preceding 
comparison continues in these figures. The gradient in the circle's amplitude is 
almost negligible in figure [231 and the negative region behind the object is again 
shallower and broader. Also in figure [2771 the dip is not as deep as in figure 1231 and 
the area behind the circle that has an amplitude close to is broader. 

The data used to reconstruct figures [2781 and [2791 encompasses 16 fold the length 
in each direction with respect to the reconstructed image, i. e. < r < 4096 
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and —2048 < x < 2048. The gradient in the amplitude of the object in figure 
I2.8l is nonexistent. An interesting difference between figure [2781 and the preceding 
images of that kind can be seen in the negative region behind the circle. The 
area has broadened, but at the right edge of the image there is a dip in contrast to 
the preceding constructions of this kind. This reminds of the reconstructions of 
the other type. Figure |2~9l shows the same pattern as the previous reconstructions 
using the approximation formula. It is worth mentioning that the area where the 
amplitude is close to almost stretches over the whole image. 
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The reconstructions in figures 12.101 and 12.1 II use double the amount of data 
in each direction relative to figures 12.81 and 12.91 i. e. < r < 8192 and 
—4096 < x < 4096. For this large amount of data both reconstructions look very 
much alike. This is due to the fact that the region where data is missing is far 
away and therefore the data from that region only has a very small influence on 
the reconstruction. 

The comparison of the preceding images suggests that using the previously de- 
rived approximative continuation of the data is superior to a continuation by 
since the reconstruction quality close to the flight track is much better. As is well 
known from experience [18], the reconstruction quality away from the flight track 
is plagued by artifacts and therefore unreliable anyway. Therefore the method of 
approximate continuation seems to be preferable. In addition the adverse effects 
away from the flight track could probably be alleviated by using an exponential 
decay as a kind of mollifier in the approximation algorithm. 
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Chapter 3 

Ghosts due to limited data 



Depending on the reconstruction algorithm there are various artifacts that appear 
in images obtained by the inversion of the spherical Radon transform, as seen in 
chapter [21 Since theorem 12.1.31 states that R is one-to-one, there should not be 
any artifacts in these reconstructions. However theorem [2.1.31 requires the data 
to be known for all x and all r > 0. As this is not the case, uniqueness is lost. 
Therefore it seems obvious that one of the main origins of these artifacts is the 
impossibility of measuring infinitely far. This phenomenon is further examined 
in this chapter. It will be shown that there exist functions with support inside the 
measured region that do not have any effect on the measured data. Functions like 
these are commonly called ghosts. A complete descriptions of these ghosts will 
be given and some examples will be shown. 



3.1 The problem of limited data 

The reconstruction formulas in theorem 12. 1 .31 corollary 12. 1 .5 1 and theorem 12.3.61 
are exact within the respective assumptions. Therefore an exact reconstruction 
without artifacts should be possible. But for each point which is to be recon- 
structed all the reconstruction formulas require data from arbitrary long distances 
with arbitrary large radii. Of course this is impossible in reality. Even in computer 
simulations this can not be achieved. 

Problems like this are known from various inverse problems. There are two kinds 
of difficulties associated with the implementation of the inversion formula for the 
spherical Radon transform. First, the data is only gathered in a discretized man- 
ner and can only be handled in a discretized manner. And secondly, the data can 
only be supported on a compact interval in opposition to the analytical model. For 
both problems there is an analogy to computerized tomography. It is known for 
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computerized tomography that the artifacts arising from the discretization have 
high frequencies and do not impair reconstruction quality if properly handled 
EOl . ETTl . Therefore this problem is not considered here, as it could probably 
be amended in a similar way and just as in the case of computerized tomography 
would probably not be as severe as the problem of limited data. The problems that 
are discussed in this chapter emerge solely from the fact that data is only collected 
over compact intervals in plane positions and circle radii. Similar problems in 
computerized tomography are the limited angle problem that is known to produce 
severe artifacts Il22l and the exterior problem, respectively. There are differences 
however, as the exterior problem still preserves the uniqueness in computerized 
tomography, although it causes some instability [23]. It will be shown in the fol- 
lowing that both limitations cause artifacts in the reconstruction of SAR-data. 
Only the two dimensional case is considered here, but the results should be easily 
transferable to higher dimensions. 

For a decent analysis of the effects of limited data it is crucial to find a set of 
orthogonal functions that represents the whole data space outside the measurable 
region because such a set allows to restrict the analysis of the effects to a well 
known set of functions. The problem is to find such a set that can also be analyti- 
cally inverted by an inversion formula for the spherical Radon transform. To this 
end the following steps are necessary. 

First, a set of orthogonal functions that are supported outside of the measurable 
region is constructed and some important properties are listed. It is shown that 
for an arbitrary, hypothetical measurement that extends over the whole half plane 
the information gathered by the projections of the data onto the set of orthogonal 
functions is sufficient to regain the data outside of some compact set that can real- 
istically be measured. Then the orthogonal functions, which span the whole space 
of data functions that could be measured outside this compact set, are inverted 
using the inversion formula in theorem 12. 1 .31 Finally some numerical examples 
are shown. 

These examples show that the problem of limited data should not be neglected. 
Since the artifacts that can arise are severe and depend on the chosen reconstruc- 
tion formula, as seen in the preceding chapter, it should be tried to find means to 
minimize these artifacts. 

3.2 Orthogonal functions and their transforms 

In the following, the maximal radius up to which the data is measured will be 
denoted with R and the start- and endpoints of the flight track will be denoted 
with — L and L respectively. Therefore the data function g(x,r) is only known 
for < r < R and —L<x<L. Now some functions and their transforms are 
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compiled. This information is later used to construct a set of orthogonal functions 
supported in R x [0, oo) \ [— L, L] x [0, R]. 

3.2.1 Definition 

J denotes the Bessel function of order zero, which is defined as 

:= D- 1 ) fc 2 2 fcfc!r(fc + 1) 

for \ arg z\ < it. 

3.2.2 Definition 

For f : [0, oo) — > R the Hankel transform of f is defined as 

oo 

J f(r)rJ (rp) dr. 
n 

3.2.3 Lemma 

Let b, b' > 0. Then 

oo 

2 f 

— / cos(frr) cos(&V) dr — bib — b'). 

7T J 


Proof: 

This is easily verified using the properties of the Fourier transform Il24ll . 

□ 

3.2.4 Theorem 

Iff : [0, oo) — > R and y/rf(r) is piecewise continuous and absolutely integrable, 
then 

oo 

f(p) = J f(r)rJ (rp) dr 
o 

exists and 

oo 

/ f(p)pJo{pr) dp = lim -(f(r + h) + f(r- h)). 
o 

Proof: 

11241 5-3, Theorem 1] 

□ 
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3.2.5 Lemma 

Let a > and 



/(*) 



J {a\/x 2 — L 2 ) for |a;| > L 
for Ixl < f . 



Then 



Q-Ly'a 2 -^ 2 



cos(x£)f(x) dx = < 



Ve-* 2 



for < |£| < a 
for |C| > a. 



Proof: 

(33 I, §17, p. 78] 

3.2.6 Lemma 

Let 



f(r) 



cos(bVr 2 —R 2 ) „ „ 

forr < R. 



Then 



cos(Ri/p 2 -6 2 ) 



rJ (rp)f(r) dr 



for p > b 



for p < b. 



Proof: 

I, 1.2, 2.56, p. 12] 

3.2.7 Lemma 

Let 



f(r) 



cos ( by/R 2 — r 2 ^j 
VR 2 -r 2 



for r < R 



forr > R. 



Then 



rJ (rp)f(r) dr 



sin (r^p 2 + b 2 



y/p^TV 
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Proof: 

$M I, 1.2, 2.55, p. 12] 

□ 

3.2.8 Lemma 

Letr, r > 0. Then the Hankel transform ofr J (r o -) is 



oo 



pJ (rp)r J (r p) dp = 5{r - r ). 



n 



The ■ is a placeholder for the variable used in the transform. 
Proof: 

This is easily verified using theorem [323] 

□ 

Now two sets of distributions will be given, {o^ nge \a £ R, b > 0} and 

{°evem °odd\ a > 0, / £ N } and it is shown that they are orthogonal. The first set 
consists of distributions that are supported only in r > R. These distributions thus 
represent the missing information due to the limitation that the reflected waves can 
only be received up to some distance R. The second set consists of functions that 
are supported only in r < R and |a;| > L. These functions therefore represent 
the missing information due to the fact that the plane travels only a limited dis- 
tance from —LioL. Note however that in this set only the information deficit is 
contained that is in addition to the first case. Therefore the constraint r < R is 
added. 

3.2.9 Definition 

LetR,L > 0. 

1. LetaeR,b> 0. Define 



2. Let a > 0,1 £ N . Define 



, . cos ( byJr' 2 —R? ) 

ab , N a(x-a) \ r ,_ Rl 



forr > R 
otherwise. 



J (ay/x 2 - L 2 ) C ° S ^J^i — for < r < R and \x\ > L 



otherwise. 
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3. Let a > 0, I G N . Define 



°odd( X ^) 



xJ (ay/x 2 - L 2 ) l ^ 2 R _ r . 2 — for < r < R and \x\ > L 
otherwise. 



3.2.10 Definition 

Define Neumann 's number e n [27\l : 

e := 1, e n := 2, n = 1,2,3, .... 



3.2.11 Proposition 

Let i?, L > 0. 



L Let o, a' G 1R, &, 5' > 0, a ^ ct', 6 7^ 6'. Then the distributions o^ n g e {x , r) 
Or'dng e (x, r ) are orthogonal to each other with respect to the scalar prod- 

oo 

uct< > 0ranae = J J f(x, r)g(x, r)ry/r 2 - R 2 dx dr. 

R R 

2. Let a, a' > 0, Z, /' G N , a ^ o', I ^ Z'. Then the functions o^ en (x, r) and 
evevk x i r ) 3X0 orthogonal to each other with respect to the scalar product 

R oo 

< f,9 >o even = J J f(x,r)g(x,r)ry/R 2 - r 2 xdxdr. 

L 

3. Let a, a' > 0, I, I' G N , a ^ a! ,1 ^ I'. Then the functions o®j d (x, r) and 
°odd ( x ' r ) 3X0 orthogonal to each other with respect to the scalar product 

R oo 

< /> 9 >o odd = J J f(x, r)g(x, r)ry/R 2 - r 2 \ dx dr. 

L 



Proof: 

Let R,L>0. 
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1. Let a, a' e E, b, b' > 0, a 7^ a', 6 7^ 6'. Then 



< a,6 a',6' > 

range' "range o ra nge 



^a,b 




OrangeK x , r )°raige( x , r)rVr 2 - R 2 dx dr 



R R 

00 




R R 



cos (by/r 2 — R? ) 

8(a — x) : — - 

V ' Vr 2 - R 2 





V 7 Vr 2 - i? 2 

00 

(<?o * y J cos (&Vr 2 - i? 2 ) cos (VVr 2 - i? 2 ) - ^ ^ 2 dx dr. 



H R 



The substitution r' = y/r 2 — R 2 and lemma [3T2.3I yield 



< <i ff e,<ie ><w 9e = (*«0(a) / cos(6r') cos(fo'r') dr' 



^(a-a')«J(6-6'). 



2. Let a, a' > 0, Z, Z' G N , a 7^ a', / 7^ /'. Then 



< o a,/ o a '''' > 

et?en' even o even 



R 00 

'evenV^i ' )" evenV^ i 




olLr, (x, r)o a '?' (x, r)r\/R 2 - r 2 x dx dr 



L 

R 00 




L 



( / \ cos (l^VR 2 ~r 2 ) / , , 

J ( aVx 2 -L 2 ) V R ^J ( a'Vx 2 - L 2 

V J Vi? 2 - r 2 V 



cos (Z'fV^-r 2 ) — 

f r v R 2 — r 2 x dx dr 

VR 2 - r 2 



R 00 



y J J (aVx 2 - L 2 ) cos (z^Vi? 2 - r 2 ) J (a'Vi 2 - L 2 ) 

cos ( Z' — V R 2 — r 2 ) = x dx dr. 

\ R J VR 2 - r 2 



L 

x 
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The substitutions r' = V R 2 — r 2 and x' = \J x 2 — L 2 and lemma 13.2.81 
yield 



R oo 



< o a ' 1 o a '' V > 

^ event ^even ^ o evcn 



f f TV 71 

I I J (ax ) cos(7— r')J (aV) cos(Z — r')x' dx dr' 





= ~5(a- a!)—5 lv 
a e { 

with the Neumann's number e/. 
3. Let a, a' > 0, /, I' G N , a ^ a', Z 7^ I'. Then 

a,i a',2' 
< °odd' °odd > o D dd 
_R 00 



J J °odd( x ' r ) 1>dd r)r v / i? 2 - r 2 ~ cir 



L 
ii 00 




xJ (aVx^I 2 ) cos ( l R VR2 r % h (a'Vx^T 2 
V J y r W~^ U V 



L 

X 

ii 00 



cos (l'^VR 2 -r 2 ) r- -1 

— rv R 2 — r 2 — dx dr 

VR 2 - r 2 x 



J J h [aVx 2 - L 2 ) cos (l^R 2 - r 2 ) J (aVx 2 - L 2 ) 

L 

cos ( l' — VR 2 — r 2 ) = xdxdr. 

V R J ^R 2 - r 2 



L 

X 



The substitutions r' = V ' R? — r 2 and x' = \J x 2 — L 2 and lemma 13.2.81 
yield 



R 00 



o,i a',;' 

< °odd' °odd > o dd 



f f 7T 7T 

/ / Jo(ax ) cos(Z— r )Jo(a X ) cos(l'—r')x' dx' dr' 




-<5(a - a')—5 w 
a ei 



with the Neumann's number e/. 



□ 
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3.3 Projection of immeasurable data onto orthogo- 
nal functions 



In the following the projections of the data onto the orthogonal function sets are 
defined, and it is shown that the data for r > R and x > L can be recovered from 
these projections. 

3.3.1 Lemma 

Let f G y and g = Rf. Then g(x, r) = g(x, —r). 
Proof: 

This is easily verified using the definition of Rf(x, r). 

□ 

3.3.2 Definition 

Let f G y and g = Rf. Define 



I. 



2 

Grange(a,b) := - < o^ nge (x, r) , g{x, r) > 0ran 



G Len( a ) ■= 7y^ a < °even{ X i r ) ,9 (x, T ) + Q {~X, r) > 0eve 



G l odd {a) := < o a J d {x, r), g(x, r) - g(-x, r) > 0odd 
3.3.3 Definition 

Let a,b G EU {— oo, oo}, a < b, and x G R. Then define the characteristic 
function of (a, b) 

1 for x G (a, b) 



X(a,b){X) 



otherwise. 



3.3.4 Theorem 

1- G range is well defined if f G 5? and g = Rf. Moreover, for r > R 



oo 

g(x, r) = J J o^'angeix, r)G range (a, b) da db. 



R 
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2. G l even (a) is well defined iff&J? and g = Rf. Moreover, if g is even in x, 
< r < R, and x > L, then 



oo °? 

g (x,r) = Y^ / o a J en (x,r)G l even (a)da. 
1=0 { 



3. G l odd (a) is well defined if f G 5? and g = Rf. Moreover, if g is odd in x, 
< r < R, and x > L, then 



oo 

oo „ 

g(x,r) = J2 / o a j d (x,r)G l odd (a)da. 

i=o i 



Proof: 

Let f E y and g = Rf. 



1. (a) First the well-definedness of G range will be shown. 



7T 

l^ 1 ]- | G range (a, b)\ 



IS 1 ] \< o a ra b nge (x,r),g(x,r) > 0r 



{S 1 
{S 1 




cos (b\/r 2 — R 2 ) /— - 

5(x — a) -^=^= g(x, r)ry r 2 — R 2 dx dr 



R R 



Vr 2 - R 2 

R R 

oo 

J rcos (b\/r 2 — R 2 ^j g(a, r) drdx 



s 1 
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With ( ri ) = y= M 

\rr}J & \y2' 



1 ^ 

\S \ — \G range (a, b)\ 



J cos (by/ y 2 - R 2 ^ f(a + y x , y 2 ) dy 



\y\\>R 

< / \f(a + yi,y2)\dy. 



\\y\\>R 



The last integral is finite, because / E 5? . 
(b) Now it will be shown that g can be recovered for r > R. 




°range( x i r)G range {a, b) da db 



R 




R 

oo 



cos (by/r 2 — R? ) 

G ran ge{a,b)X(R,oo){r)5(x - a) da db 

Vr — R z 




7T 




X(R,oo)(r) - / / S(a- x')r' cos ( bV r' 2 - R 2 



o R 



R R 



cos (by/r 2 — R? ) 
xg(x,r)drdx \5(x — a) , c?ac?6 



oo / oo 



Vr 2 - i? 2 

,oo)( r ) J ~ \ J r ' cos r' 2 — R 2 ^ g(x, r') dr' 



o \R 



cos (by/r 2 — R? ) 

x y 1 db 

y/r 2 - R 2 



With the substitution r" = y/r' 2 - R? and [g (x, y/- 2 + R 2 )] indi- 
cating the cosine transform of r"g (x, yf r" 2 + R 2 ) with respect to r" 
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it follows that 



oo 

j j °range( X i r ) G range(a,b) da db 



R 



OO / OO 

X(R,oo)(f) J ~ij cos (br")r"g (x, V r" 2 + R? j dr" 
o \o 



cos {b\Jr 2 — R? ) 

x - — '- 

Vr 2 - R 2 

X(R,oo)(r) 

Vr 2 — R 2 J 
o 



db 



g[x,yp*T&\ (6) 



cos ^feV r 2 — R? j d& 



X(B,oo)( 



==Vr 2 ~ R 2 9 (x, \l Vr 2 - R 2 * + R^j 



X(R,oo)(r)g (x,Vr 2 - R 2 + R 2 j = X(R,oo)(r)g(x,r). 



2. Assume without loss of generality that g is even in x. 



(a) First the well-definedness of G even will be shown. 



eia 



\ Sl \\\< ° a eL( x i r),g{x, r) + g(-x, r) > c 



IS 1 ] 



R oo 



J ( aV x 2 — L 2 



o L 



cos (l^VR 2 ~r 2 ) 
VR 2 - r 2 



x (g(x, r) + g(—x, r))rv R? — r 2 x dx dr 
\S l 



J xJ (aV x 2 — L 2 

L 
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R 



J xJ (a^x 2 - L 2 ) J cos (^V# 2 - r 2 



x (//(x + re.n,)^,,)) *d. 



<cj j jMb + ^teistf,,)*- 

S 1 L 

with an appropriate constant C > 0. The change in the order of inte- 
gration is valid and the last integral is finite because / G 5? . 

(b) Now it will be shown that the even part of g can be recovered for 
< r < R and Id > L. 



Yl / ° a eLn(x,r)G l even (a)da 

OO °° 

= (x(-oo,-L)(x) + X(L,oo)(x)) X(o,R)(r) Y / J ° (aVx 2 - L 2 ) 



X 



X 



cos (l^VW^) / €l f , ( f ^— r2 

w - r 2 ^ / ax j ° - 



J C ° S Oi?^ ^ ~ r ' 2 ) r '^( x '' d(~ x ' ' r 0) ^ J ^ a 
o / 

OO 

OO „ 

(X(-oo,-z,)(z) + X(l,oo)0)) X(o,H)(0 ^2 ho (aVx 2 - L 2 

1=0 { 

cosQl^g) / e, ? , 
V# 2 - r 2 \ i?7 V 



X 



x 



t\/ x' 2 — L 2 



R 



J C ° S 0^2^ ^ ~~ r ' 2 ) r '9( x ' ,r '^ dr' dx' J da. 



The substitutions x" = V x' 2 — L 2 and r" = \/ R 2 — r' 2 yield with 
gi iys/x" 2 + L 2 ) representing the Fourier coefficients corresponding to 
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r" g (V ' x" 2 + L 2 , \/R 2 — r" 2 ) and (V - 2 + L 2 ) indicating the Hankel 
transform of gi {\/ x" 2 + L 2 ) with respect to x" 

oo °° 

J] / o^ en (x,r)G^ en (a)rfa 
1=0 { 

oo ~ 

= (x(-oo,-l)(x) + X(L,oo)(x)) X(o,R)(r) ^2 J J ° (aVx 2 -L 2 ) 







cos (lj,VR 2 -r 2 ) e t f . „. 
V f — ~tt / ax"J ax" 



X 



/ C ° S OtT ") ( v/x " 2 + L2 ' Vi? 2 - r" 2 ) dr"dx"da 



o 



X 



oo oo 

y aJo (aV x 2 — L 2 j / x"Jo(ax")gi (^V x" 2 + L 2 j dx" da 





(x(-oo + X(L,oo)(a;)) X(o,fl)(r) 



cos (^I—VR 



X 



aJ (Wx 2 - L 2 ) # (V- 2 + L 2 ) ( a) c?a 
o 

(X(-oo,-L)(x) + X(L,oo)(g)) X(o,-R)(y) 

v^ 2 - r 2 

x ^cos (l^U/R 2 -r 2 ^ gi N Vx 2 - L 2 " 2 + L 2 ) 

(X(-oo,-L)(x) + X(L,oo)(x)) X(0,R)( r ) 



r 



Vr 2 - 

x ^cos {l^\/R 2 — r 2 J aj(x) 
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(X(-oo-L)(x) + X(L,oo)(x)) X(o,R)(r) 

VR 2 - r 2 



x V R 2 — r 2 g l x, 



R 2 -VW 



= (x(-oo,-l)(x) + X(l,oo)(x)) X(o,B)(r)g(x,r). 
3. Assume without loss of generality that g is odd in x. 
(a) First the well-definedness of G Q dd will be shown. 

eia 



\ S \ < o a j d {x,r),g{x,r) - g(-x,r) 



°odd 



1^1 



R OO 

/• /• / -\ cos (l^VR 2 -r 2 ) 

L 



x r) — g(—x, r))r\/R 2 — r 2 — dx dr 

x 

OO 

IS 1 ! y a;Jo (a^x 2 - L 2 ) 

L 

x / cos (Z — \/ R 2 — r 2 ) —rg(x, r) dr dx 
J V R J x 

o 

oo R 

J J (Wx 2 - L 2 ) y cos (^v 7 ^ 



2 2 I y» 



R oo 



< 



C j j J \f(x + r£,rr))\dx dSx(£,r))dr 



S 1 L 



with an appropriate constant C > 0. The change in the order of inte- 
gration is valid and the last integral is finite because / E y. 

(b) Now it will be shown that the odd part of g can be recovered for < 
r < R and \x\ > L. 
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CO ~ 

Yl / odd(^ r ) Gl odd( a ) da 
1=0 { 

oo 

oo „ 

= (x(-oo-l)(x) + X(l,oo)(x)) X(o,R)(r) Y / xJ ° ( a ^ x2 ~ L2 

i n v 







x f — / ax J ayx' 2 - L 2 

VR 2 - r 2 \ 2RJ l V 

R \ 

x y cos (l^\/R 2 — r /2 j r'-^j(g(x', r') — g(—x', r')) dr' dx' J da 
o / 

(X(~oo-L)(x) + X(L,oc)(x)) X(o,R)(r)^2 / xJ ° (aVx 2 - L 2 







X 



/ oo 

v^ 2 - r 2 \rj v y 



X 



o 



With the substitutions x" = y/x' 2 — L 2 , r" = R 2 — r' 2 , 
gi (v 'x" 2 + L 2 ) indicating the Fourier coefficients corresponding to 

r" g \\/x" 2 + L 2 , y/R 2 — r" 2 ) and i^, L % g\ (V - 2 + L 2 ) representing 

the Hankel transform of ^ X ,} +L i 9i (Vx" 2 + L 2 ) it follows that 

oo ~ 

^ / °odd( X i r ) Gl odd( a ) da 
1=0 { 

oo 

oo „ 

= (x(-oo,-L)(a?) + X(L,oo)(^)) X(o,ji)0*) X] / xJ ° (aVa; 2 - L 2 

7 n J 



cos {l^y/W^r 2 ) e t 



i=o 

oo 

a Jo (ax") 



- r 2 i?7 Vx" + L 2 



o 
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R 

X 



J cos (^i? r ") r " ^ x " 2 ^ ^ 2 ~ r " 2 ) 



00 cos (liLyjR? — r 2 ) 
(X(-oo,-l){x) + X(L,oo)(a;)) X(o,ii)(r) V f 

00 00 
y a Jo (a\/ x 2 — L 2 ^ ^ x"Jo(ax?') 


=a> fV x" 2 + L 2 ^ dx" da 

Vx" + L 2y V y 

(X{-oc,-L)(x) + X(L,oo)(x)) X(0,R)( r ) fp 2 _ 2 



X X 



X 



X X 



J aJ (oVx 2 - L 2 j — = = == g ; (V- 2 + L 2 ) (a) da 



fa-oo,-*) (*) + X(L,oo) (*)) (r) x £ cog (jl.JW=?) 



x ^=p; V v^ 2 - L 2 + L 2 



Vx 2 - L 2 + L 2 

(x(-oo,-L)(x) + X(l,oc)0)) X(o,fl)(r) 



VR 2 - r 2 

00 

X 



00 

^^cos V i? 2 — r 2 J <ft(x) 



(X(-oo,-L)(g) + X(L,oo)Qe)) X(0,R)( r ) 

Vi? 2 - r 2 
x V/2 2 - r 2 fx, 

(x(-oc,-L) (x) + X(l,oo) (x)) X(o,R) (r)g(x, r) . 



□ 

Note that this theorem uses essentially the cosine and the Hankel transform and 
it looks like an overly complicated formulation. However it is essential to choose 
the functions like this to compute the inversions analytically. 
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3.4 Reconstruction of the orthogonal functions 

The previous chapter showed that the missing data from the unmeasurable region 
can be projected onto the set of orthogonal functions {o^ nge : a G R, b > 0} U 

{Og4n> °o'dd : a > 0, / G N } and can be retrieved again. In the following it 
is therefore sufficient to consider only the orthogonal functions to examine the 
error that arises from limited data. Now the reconstructions of these functions are 
performed. 

3.4.1 Definition 

Letf G y(R n x R), a G R", and b G R. Then (r {a>b )f)(x,y) = f(x + a,y + b). 

3.4.2 Theorem 

1. Let a G R, b > 0, and 



g(x,r) = o\ 



a,b 



range 



Then Rf = g with 




and 




where H y refers to the Hilbert transform in y. 



2. Let a > 0, / G N , and 



g(x,r) = o a J en (x,r). 
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Then Rf = g with 



/«„) * 



' g-Ly/a 2 -£ 2 



x < 



sinfl/-^/^ — 



/or < |f | < a 



/or |f | > a 



and 



2 3 



it y dy 



R 



»(l%Vfl 2 -(s-t) a -y a ) 



/ (oVt 2 - £ 2 ) /or |t|>L 
/or |t| < L 

\ 



/or (x - t) 2 + y 2 <R 2 



for (x - t) 2 + y 2 > R 2 J J 



eh 



3. Let a > 0, I G No, an J 
TTzen i?/ = g with 

m v) = - 



g (x,r) = o a Q 'j d (x,r). 



1 d sin (Ry/? + ri 2 + {l%) 2 ) 



for < |f | < a 



' g-L\A 2 -? 2 



x < 



sin(L-y/^ — 



/or |f | > a 



and 



2 9 

-XHy — 



7 (a Vt 2 - £ 2 ) far \t\ > L 



r 







for \t\ < L 



for (x - t) 2 + y 2 <R 2 



x 



sjR 2 -{x-t) 2 -y 2 

for (x - t) 2 + y 2 > R 2 
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Proof: 

1. Let a E R, b > 0, and 



g{x,r) = o^ nqe (x,r). 



range \ 

Then the three dimensional Fourier transform of g is given by 

oo 
R 

and with lemma UX6l 



V Z7T 



cos(rV^) 

1 / 2 h 2 for p>b 



for p < 6 

Therefore, according to theorem 12X3] 

cos(i?-y/$ 2 +r; 2 -b 2 ) 



for j£ 2 + 7? 2 > 6 

/ 8tt 



for + r/ 2 < 6. 



Hence 

-1 = 

2tt 



f(x,y) = ±- I I f&fM f{tri)d-qd£ 




R R 

1 



2tt 




R R 
1 



X 



cos f R^fi 2 + r] 2 -b 2 
= \v\e-^ V ^ 2 / -x (& 2 ,oo)(e 2 + ^ 2 ) | <M£ 

7T +7] 2 — ¥ 



1 --T(- a , )KyS- I 




R R 



cos (iVf 2 + ^ 2 - b2 ) 



X(fc 2 ,oo)(^ + ?7 ) dr ld£, 



+ t/ 2 - b 2 

1 9 7 1 cos ( R ^ p2 ~ 62 ) 
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With lemma 1X2 .61 and because the Hankel transform is its own inverse, this 
leads to 



8tt 



T(-a,0)Hy 







cos( by/ x 2 +y 2 —R 2 



for yjx 2 + y 2 > R 



y/x 2 +y 2 -R 2 

for yjx 2 + y 2 < R 



=H„— — 



dy 

cos (by/(x-a) 2 +y 2 -R 2 ^ 



for yjjx - a) 2 + y 2 > R 



>&F V dy 
2. Let a > 0, / € N , and 



y/(x-a) 2 +y 2 -R 2 

for ^ - a) 2 + y 2 < i?. 



g(x,r) = o a J en (x,r). 
Then the three dimensional Fourier transform of g is given by 

oo 

= / / e ~ l ^ ri o( r P)9( x i r ) drdx 



R 



and with lemmas I3.2.5l and l3.2.7l 



2 sin (Ry/pP + jl*)*) 



-L^/at-t 2 



y/a^e 
sm(Ly/(, 2 -a? 



for < |£| < a 
for |f | > a. 



Thus, with theorem [2X3] it follows that 



i H sin(iVe+^TT7IF) 



x < 



for < |£| < a 



for |£| > a. 



Therefore 



R R 
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Applying the Fourier convolution theorem leads to 

/( ^ )= (2^ H 4/ 

R 

/ / „ r./T2 t2 

for < |£| < a 



( 



R 



v^i 2 

sin (L ^ -a 2 ) 

v^ 2 



for |£| > a 



x 




• / M . sin (i? A A 2 + r? 2 + (/£) 2 ) 



R R 



ix^ y dy 



R 



' g-Z,\/a 2 -? 2 



cos(t£) < 



V^i 2 



V 



for < |£| < a 



for |£| > a 



With lemmas [3 .2.5l and r3.2.7l and bearing in mind that the Hankel transform 
and the cosine transform are their respective inverses this leads to 



f(x,y) 



dy 



r 



J (aVt 2 -L 2 ) for \t\ > L 
for \t\ < L 



x 







for - tf + y 2 <R 2 



for (x - t) +y > R J J 



dt. 



3. Let a > 0, / G N , and 

g (x,r) = o a j d (x,r). 
Then the proof is analogous to 2. using x sin(x£) 



-^cos(z£)- 



□ 
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3.5 Numerical simulations 



To imbue the unwieldy formulas from theorem 13.4.21 with life, in the following 
some of the ghosts - functions that are in the null space of the measurement oper- 
ator due to limited data - derived in this theorem are shown. All of the following 
examples can therefore be added to the reflectivity function / without changing 
the measured data. In all cases it is assumed that L = R = 1, and the region 
(x, y) G [0, 1] x [0, 1] is shown. At first, examples for g G {o^ nge : a G R, b > 0} 
are shown, then examples for g G {o a J; en : a > 0, / G N }. The examples for 
g G {o a ^ d : a > 0, 1 G N } are omitted because they differ from the examples of 
g G {Og^ en : a > 0, I G N } only by an additional factor x. 



3.5.1 Reconstructions for g e {o^f : a e R, b > 0} 

Now two series of images are shown to demonstrate the effects of the parameters 
a and b for ghosts derived from functions g G {o^ nge : a G R, b > 0}. These 
are the artifacts that can arise due to the limitation that the echoes from the emit- 
ted waves can only be received up to the limited distance R. In the first series 
b is constant at 0.25, and a is varied from —0.1 to 1.2. That range is sufficient 
because reconstructions for different values of a can be obtained via translations 
in x, as can be seen from the formula in theorem [3.4.21 1. In the second series 
a is constant at 0.6, and b takes the values of 0.25 and 1. Since the functions 
g G {orange : a G R, 6 > 0} all have a singularity for r = R that would usually 
not be reflected in real data and that would only obscure the interesting features, 
only the differences o^ nge — are reconstructed. Thereby the singularity is 

removed and that should lead to ghosts similar to those that can be expected in 
reconstructions from real data. The maximum displayed in the following images 
is capped at a reasonable level to enhance the visibility of the features. Otherwise 
a couple of large peaks would distort the image. These peaks are located at the 
points where the course of the depicted circles is close to a grid point. It should 
be noted that the following reconstructions are only approximative because the 
Hilbert transform from theorem [3.4.21 does only exist in the distributional sense, 
even using the differences mentioned above. This should not be a problem how- 
ever since for real data the integration with the kernel g r as in theorem l3.3.4l should 
lead to a regular function. 

Figures l3Tl to l3.4l each show a reconstruction of o^ n 2 g 5 e —o^ nge . Figure lBTTl displavs 
the result for a = —0.1, figure I3T21 for a = 0.3, figure [331 for a = 0.8, and figure 
13 .41 for a = 1.2. It can be clearly seen in figures [3TT1 to 13 .41 that the missing data 
in this parameter range causes a circular artifact that is centered at (a, 0). These 
artifacts remind of the artifacts seen in chapterEl so this is a possible explanation. 
It should be noted that the reconstructed functions are not for (x — a) 2 + y 2 > R 2 
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although some reconstructions could lead to this assumption, but this is only an 
effect of the overshadowing singularity. In the cross section of figure 13.41 it is 
clearly visible that the reconstruction is greater than for (x — a] 



y 2 >R 2 





Figure 3.5: Reconstructed / for g = Rf = o\ 



0.6,0.25 _ „0.6,0 
range range 



f(hv) 




Figure 3.6: Reconstructed / for g = Rf 



„0.6,1 _ „0.6,0 
u range u range 



Figures 13.51 and 13.61 each show a reconstruction of o° a 6 ^ e — o^° ge . Figure 13.51 
delineates the result for b = 0.25 and figure [3T61 for b = 1. As can be seen, apart 
from the scaling figures 13.51 and 13.61 are very similar. They depict a circle with 
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center (0.6, 0) and radius r. It can be seen in the cross section that the amplitude 
changes due to the variation in the second parameter. 



3.5.2 Reconstructions for g e {o a J; en : a > 0, / G N } 



In the following two groups of images, functions g G {o^ en : a > 0, 1 G N } 
are shown. The first group, for I = 1, comprises examples for a G {1,4, 16} 
as does the second group, however for I = 16. As in the preceding subsection, 
the functions g G {o^ en : a > 0,1 G No} all have a singularity for r = R 
that would usually not be reflected in real data and that would only obscure the 
interesting features. This problem is solved in the same way so that the shown 
examples should lead to ghosts that can be expected in real data. As in the previous 
subsection it should be noted that only an approximation is computed. But again 
this should not be a problem as the same circumstances hold. 

f(.bv),f(iv) 





1 












Figure 3.7: Reconstructed / for g = Rf = o, 



i.i 



4,0 



Figures [3771 to 13.91 each show a reconstruction of o a e £ en — o a J? en . Figure [3771 depicts 
the result for a — 1, figure [3T8l for a = 4, and figure [3791 for a = 16. Figures [3771 
to 13.91 again show circular artifacts, but of a different kind. These artifacts show 
the biggest variation in amplitude around x 2 + y 2 = L 2 . However the artifacts 
clearly extend into the region x 2 + y 2 < L 2 . Figures 13.71 and 13.81 remind of 
artifacts that appear in reconstructions of objects close to the edges of the flight 
track whereas figure 13.91 exhibits a higher frequency phenomenon that was not 
encountered so far. In comparing figures 13.71 to 13.91 one can conclude that with 
larger a the amplitude of the artifact gets smaller and the frequency of the artifacts 
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becomes larger. This is understandable as the parameter a affects the frequency in 
the Bessel function. 




Figures 13.101 to 13.121 each show a reconstruction of o a e ^ n — o a e ^ en . Figure 13.101 
displays the result for a — 1, figure [3TTT1 for a = 4, and figure I3T21 for a = 16. 
For these figures the same holds as for figures 13.71 to 13.91 Again, similar to 
the preceding subsection, an increase of the second parameter leads to a higher 
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amplitude as demonstrated by figures 13.101 to 13. 121 Additionally, the artifacts 
seem to get sharper. 

No reconstructions for the functions g e {o^j d : a > 0,1 G N } are shown since 
they would differ from the examples in this subsection only by an additional 
gradient in x direction as can be clearly seen in theorem [3A2l 

The computed ghosts partially show a close resemblance to the artifacts known 
from reconstructions. Unfortunately the analysis of this chapter shows that the 
common and unavoidable problem of limited data is not restricted to high fre- 
quency artifacts. This urges a thorough examination of how these artifacts can 
be avoided, e. g. by a regularization using information from the physical back- 
ground. The knowledge about the nature of the ghosts that was derived in this 
chapter should hopefully simplify this task. 
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Chapter 4 



A new approach to invert the 
spherical Radon transform 



With the ideas from the last chapter, it is possible to reconstruct the recoverable 
part of / using orthogonal functions. In the following, a set of two dimensional, 
orthogonal functions with a compact support will be introduced. Then it will be 
shown that for functions that fulfill the properties of a measurement the infor- 
mation gathered by the projections of the data onto the orthogonal functions is 
sufficient to regain the data. Finally, the inversions of the orthogonal functions 
will be calculated. With these results an alternative way to reconstruct the images 
is obtained: First, the data is projected onto the orthogonal functions, and then the 
inversions of the orthogonal functions are summed using the coefficients obtained 
from the projections. 

4.1 Orthogonal functions and their transforms 

Now a set of functions is defined, and it is shown that they are orthogonal. 
4.1.1 Lemma 

Let L > 0, < x < L, and a > 0. Then with the Bessel function of order zero J 



L 




1 



cos 



VL 2 - x 2 



cos 



(x£)dx=^J (Ly/tf + P). 







Proof: 

m I, §5, p. 30] 



□ 



61 



4.1.2 Definition 

Let R, L > and k, I e N . Define 

1. 

cos {k\^L 2 -x 2 ) cos (l^VR 2 ~r 2 ) 



hven( X i r ) ■= X(0,L)( X )X(0,R)( r ) 



VL 2 - X 2 y/K- 



cos (kfVL 2 -x 2 ) cos (l^VR 2 ~r 2 ) 



iodd&r) '■= X(o,L)( x )X(o,R)( r ) x - 



VL 2 - x 2 \[W^ 



4.1.3 Proposition 

Let R, L > 0, fc, fc', I, V E N , fc ^ k', and I ^ Then 

i. the functions ^en( x ; r ) iev'en( x i r ) are orthogonal to each 

other with respect to the scalar product < f,g >i euen = 



J J f(x,r)g(x,r)r\/R 2 — r 2 x \J L 2 — x 2 dxdr. 
o o 

2. the functions i^i(x,r) and iodd( x i r ) we orthogonal to each 

other with respect to the scalar product < f,g >i odd = 

R L 

J J f(x, r)g(x, r)ry/R 2 — r 2 ^ L2 ~ x2 dxdr. 



Proof: 

Let R,L>0 and k, k', I, V e N . Then 



1. 



R L 



< i even even >ie„e«— J J i>even( X i r )^even( X J r ) r VR 2 — T 2 X V L 2 — X 2 dx dr 


L cos {k^VL 2 -x 2 ) cos (l^VR 2 ~r 2 ) cos (k'l^L 2 - x 2 ) 








y/L 2 - x 2 VR 2 - r 2 VI 



2 — x 2 



x — f= - 7" y/ R? — r 2 x V L 2 — x 2 dx dr 

VR 2 - r 2 

y y cos (jc—VL 2 — x 2 ^j cos — V-R 2 — r 2 j cos ^k'—y/L 
o 

cos (/' — V R 2 — r 2 ) = . — (ir. 

V R J Vi? 2 - r 2 VL 2 - x 2 



2 — x 2 





X 
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With the substitutions r' = \/R 2 — r 2 and x' = \J L 2 — x 2 

R L 

< ievem ievL >w„ =11 cos(kyx') cos(Z-Jr') cos(k'yx') cos^'^V) dx dr' 



L v R ' L y R 







— Okk'—Ow 



with the Neumann's numbers e^, e\. 

2. 



R L 




k I -k' V I I -k I / \ -k' I' i \ /tto o — X 2 7 7 

nM^nM >;„,,= / / i^(a?, r)z J A (x, r)rVR 2 - r 2 dxdr 



<l oddi l odd > iodd~ 

Jb 



R L 




cos (kfVL 2 -x 2 ) cos (ij^VR 2 -r 2 ) cos (k'l^/L 2 -x 2 ) 







VL 2 - x 2 VR 2 - r 2 v 7 !^ 



x — R — — r\/ R 2 — r 2 — — X da; dr 

Vi? 2 - r 2 x 

R L 



J J cos (k—^/L 2 — x 2 ^ cos (I — VR 2 — r 2 ^ cos (k'—V L 2 — x 2 ^ 

o 

cos il' — y/R 2 — r 2 ] = , X =dxdr. 

V R J ^R 2 - r 2 VL 2 - x 2 



o o 

X 



With the substitutions r' = v 1 R 2 — r 2 and x' = \J L 2 — x 2 

R L 

•fcj -k'V f f - 7T „ ,.7T „ ,7T „ 



< i oddi i odd >i odd = / j co8(k-x')cos(l-r')cos(k'-x')co8(l'-r')dx'dr' 







— — Okk'—Oll' 

tk ei 

with the Neumann's numbers e^, ej. 

□ 



4.2 Projection of the data onto the orthogonal func- 
tions 

It will be shown that the measured data can be projected onto the set of orthogonal 
functions and that it is completely recoverable from these projections. 
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4.2.1 Definition 

Let f G y, g = Rf, and let e^, ei denote the Neumann's numbers. 



GeL ■■= t|t| < i k J m (x,r),g(x,r) + g(-x,r) > ievm 



G oL ■= 777% < *tdd( x > r )>9(x,r) -g(-x,r) > iodd 



2LR 



4.2.2 Theorem 

1. (jgien * 5 we ^ defined, and if f G , g = Rf, < r < R, < x < L, and 
g is even in x, then 



g(x,r) ^ ] ^eveni. X y r ) G even- 



k,l=0 



2. G 'L is well defined, and if f G SP, g = Rf, < r < R, < x < L, and g 



odd 

is odd in x, then 



oo 

g( X ,r) = J2 ik odd( X ^) G oL- 

k,l=0 



Proof: 

Let / G y, g = Rf, < r < R, and < x < L. 



1 . G^ en is obviously well defined since the integrals are finite. Assume with- 



64 



out loss of generality that g is even in x: 

oo 



'even 

fcj=0 



~ cos (kf VL 2 - x 2 ) cos (IWR 2 - r 2 ) 

WOW*)L VL2 _ x2 VR2 _ r2 

x |jr| < ieLn(x,r),g{x,r) + g(-x,r) > ieven 
* / / C0S (^v^^) cos 

^ 

x x'r'(g(x' , r') + g(—x', r')) dx' dr' 



X(o,^)(r)x(o,,)(x)^ Vl2 _^ 2 



£fc£[ 
LR 



With the substitutions x" = VL 2 — x' 2 and r" = V R 2 — r' 2 
and g k> i representing the Fourier coefficients corresponding to 
f%x"r"g(VL 2 - x" 2 , VR 2 - r" 2 ) it follows that 

oo 

Ei k ' 1 (r r)G k ' 1 
''evenV^i ' ) lu even 

X(oM r )X(o,L)( x ) 2^ T7W=, 72 T7WTT2 



k,l=0 



Vr 2 



L R 



x I -j^ / / cos \k^-x") cos (l-^-r" 



LRJ J \ L J \ R 

o o 



65 



x x"r"g (^V L 2 — x" 2 , V R 2 — r" 2 ^ dx" dr" 



VL 2 - x 2 VR 2 ~ r 2 

oo 

: co 

k,l=0 

1 1 



cos ^A;— V L 2 — x 2 j cos V i? 2 — r 2 

fc,i=0 

X(o,^)(r)x( ,i)(^ 



\/£ 2 - a; 2 V^ 2 - r 2 
x VL 2 - xVi? 2 - r 2 g U/ ' L 2 - VU^ 2 , \/ R 2 - VR 2 - r 2 2 

X(o,R)(r)x(o,L)(x)g(x,r). 



2. G^' dd is obviously well defined since the integrals are finite. Assume without 
loss of generality that g is odd in x: 



k,l 
odd 



k,l=0 



v Wv M V /os^fyT^jcos^x/i^P) 



k,l=0 



x < r),g(x, r) - #(-x, r) > iodd 

(r)X(o,,) (x) X, x - ? =_ 7 === 

^//cos(^Vra)cos(^ 



x r'(g(x', r') — g(—x', r')) dx' dr' 



2 _ r /2 
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X(0M r )X(0,L)( X ) 2^ X 2 ./P2— 3 



VL 2 - x 2 VR 2 - 
J J cos (k^VL 2 - x' 2 ) cos (^Vi? 2 - r' 2 ) 



fc,i=0 

L i? 

Li? 





x r'g(x', r') dx' dr' 



2 _ r i2 



With the substitutions x" = \/L 2 — x' 2 and r" = y/R 
and ^ fc / representing the Fourier coefficients corresponding to 



LR yJU^i 



"g(VL 2 - x" 2 , VR 2 - r" 2 ) it follows that 



E-k,l / \s-ik,l 
l odd\ X i r ) U odd 

k,l=0 



v frW M V r ^HVU^) cos (i^VW^) 

x(oM r )x(o,L)( x ) 2^ x — TF^P — 



X 



k,l=0 

L R 

k 1 J J cos {^Y^'^j cos (^ r ") 



Li? 



x - X r^fVL 2 -x" 2 , VR 2 - r" 2 ) dx" dr" 
VL 2 - x" 2 V 

X{o,K)(r)x(p,L)(x)x 





II 



VL 2 - x 2 VR 2 - r 2 



cos (kjVL 2 - x 2 ) cos (^v 7 ^ 2 

fc,/=0 

X(o,i?)(^)X(o,L)(a;)a;- 



k,l=0 

1 1 



a/L 2 - X 2 VR 2 



VL 



2 — X 2 



=\/R 2 — r 2 



L 2 -^L 



x g(yj L 2 - VL 2 - x 2 , \JR 2 - VR 2 - r 2 2) 
X(o,fl)(r)X(o,£)(a;)^(ar,r). 



□ 



67 



4.3 Reconstruction of the orthogonal functions 

As seen above, the measurable data can be projected onto the set of orthogonal 
functions {ieven^tdd '■ k,l £ N } and can be recovered again. In the follow- 
ing it is therefore sufficient to perform the reconstruction only for the orthogonal 
functions. 

4.3.1 Theorem 

1. Let k,l G N and 

g{x,r) = i^ en {x,r). 

Then 

sin (iiy(Zi) 2 + £ 2 + 77 2 ) 



and 



I vv-* f° r \t\< L 

\ y otherwise 



x 



costey/V-fr-V-V 2 ) , , ^2,2 D 2 1 \ 

— v , = = ' for (x - t) 2 + y 2 < R 2 



dt 



otherwise J J 

and H y refers again to the Hilbert transform in y. 
2. Let k, I G No and 



Then 



f&v) = -\h\v\^Jo[Lj(k-) 2 + Z 



and 

/ ( cos(fcf Vi 2 -* 2 ) 

/(z,2/) 





OlV^-M 2 -^ 2 ) , / + N2 , 2 D 2 1 \ 

v / ' for (x - tY + y 2 < R 2 

otherwise J J 



dt. 



68 



Proof: 



1. Let k, I e IN and 

g(x,r) = i h J en (x,r). 
Then the three dimensional Fourier transform of g is given by 

oo 

9(Z,p) = ^7= J J e~ txi rJ (rp)g(x,r) dr dx 



JR. 

and with lemmas l4.1.1l and [3.2.7l 



7T. 2 \ Sin (Ry/{l%)*+f?) 



g(^p) = J- J Q LJ(k-) 2 + e , 

With theorem f2. 1.31 this results in 



T Hh[LJ{kly+e 



sin (R^/(l^ + e + V 2 ) 



Therefore 

JR. JR 

and applying the Fourier convolution theorem leads to 



R \R 

R R 



d f I f , ^ ( , 7T 



2%^/ l/cos^OJo^Cfc^ + e ) # 

R \0 



■DC 



v sin (RJ(l^) 2 + p 2 ) 
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With lemmas 14.1.11 and 13.2.71 and because the Hankel transform and the 
cosine transform are their respective inverses, this can be written as: 



f(x,y) 



2 y dy 



( ( cos(fc fVL 2 -i 2 ) 
VL 2 -t 2 



TR 



X 



V / J? 2_ (a ,_ t) 2_ y 2 




for (x - t) 2 + y 2 <R 2 



dt. 



otherwise J J 



2. Let fcjeNn and 



9{x,r) = i k J d {x,r). 



Then the proof is analogous to 1. using x sm(x^) = — J| cos(x£). 



□ 



4.3.2 Remark 

Geven m & G odd C 

transform as described in f28l/ . As according to theorem \4.2.2\ 



GeLn m ^ G k odd can be computed via a fast two dimensional non-equidistant cosine 



fc,Z=0 



if follows with the results of theorem [4.3. 1\ that 



fc,2=0 



j \/L 2 -t 2 

Vio 



for k| < L 



otherwise 



cos (l^B?-{x-t) 2 -y 2 ) ^2,2 D 2 1 \ 

— ^ for (x - t) 2 + y 2 < R 2 



y/R 2 -(x-t) 2 -y 2 







dt. 



otherwise J ; 



This summation of the precomputed solutions of the orthogonal functions with the 
appropriate weights G k J; en and G odd is not faster than the direct implementation of 
the inversion formulas presented so far. However two improvements are possible. 
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1. This inversion formula is adapted to the fact that the measured data is lim- 
ited. Therefore it is probably easier to identify means to diminish the arti- 
facts caused by this limitation. 

2. Depending on the application, the coefficients representing higher frequen- 
cies are probably prone to noise. Therefore it could be worthwhile to neglect 
them in the summation to obtain a faster algorithm since these coefficients 
do not contain much dependable information. 
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Chapter 5 

Tackling the left-right-ambiguity 



The so called problem of the left-right ambiguity is one that exists since the in- 
vention of SAR and continues to be an issue, albeit alleviated, till this day ||29l . It 
refers to the fact that because of geometrical reasons the reconstruction formulas 
in this field are not able to distinguish between points that lie symmetrically to the 
left and right of the flight path. At first this problem seemed inevitable lfT3l . lfT2ll . 
but later, approaches to solve it were proposed. The method of beamforming (SJ, 
for example, tries to overcome the ambiguity by directing the energy emitted from 
the antenna as much toward one side of the flight track as possible. In this way 
the echo from the illuminated side is much stronger than the echo from the other 
side. Unfortunately these weak signals still lead to shadows in the reconstruction 
for highly reflecting objects. A different solution was proposed in ll30ll . The idea 
consisted of using data from two parallel flight tracks. However, singularities aris- 
ing from noisy measurements present a serious problem in this approach. 
Inspired by the idea of gathering and then combining data from two slightly dif- 
ferent positions, two post-processing formulas are given in the following that ma- 
nipulate the data measured by an airplane equipped with two or more antennas. 
It is shown that it is possible to recover the odd part of the reflectivity function 
from measurements with at least two antennas. The reconstruction of the origi- 
nal image, including the odd part, is difficult for only two antennas, but it will be 
shown that since the low frequencies are attenuated in the inversions of the spher- 
ical Radon transform, it is possible. As will be seen, for a good reconstruction 
quality it is nevertheless recommendable to use more than two antennas. This has 
the additional advantage that no regularization is needed and therefore an exact 
formula can be used. 

The chapter is structured as follows. The first section will present analytic consid- 
erations, which prove that it is possible to retrieve the original, asymmetric image 
given at least two even images. This corresponds to a measurement with at least 
two antennas. In the next section the two post-processing formulas are introduced. 
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The first requires a sophisticated regularization whereas the second needs at least 
three source-images as input. Afterwards the results of numerical simulations are 
displayed. 

5.1 Using two data sets 

At first some definitions are necessary. Then the important trick is shown sepa- 
rately, which makes the following central theorem fairly easy. At the end some 
conclusions derived from the theorem are mentioned. 

5.1.1 Definition 

Let f e y(R n x R) and b e R. Then 

1. f F ' J and f I,F denote the Fourier transform of f in the first n variables and 
the last variable respectively, f 1,1 = f and f F,F = f. 

2. / (c ' /} (x, y) = f{-x, y) and / (/ ' c) (a;, y) = f{x, -y). 

3. f§{x,y) = l\f(x,y + b) + f(x,-y + b)] 

5.1.2 Lemma 

Letf e y(R n x R) and be R. Then (r (0 , b) f) {I ' c \x, y) = rp,,^)/^. 
Proof: 

{T {0 , b) f) {I > c \x,y) = (T (0)6 )/)(x,-y) = f(x,-y + b) 

= fV' c \x,y-b)=T {0 ,_ b) fV> c \x,y). 

□ 

Inspired by the ideas in QUI , the following theorem is formulated. 

5.1.3 Theorem 

Let f e y(R n x R), 6, rj E R, and x 6 R n . Then 

sHbri)f> F {x,rj) = l -\r b - r i0 ^ b) f^ F (x, V ). 

Proof: 

fmibrfif'Sfarj) = ^f'^x^-^f'^x,^} 
= Trrfao, b)f - r (0 , - b) f} I>F (x, rj) 
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^[W - r (0 ,- 6) / + (r (M )/) (/ ' c) - ( W) (/ ' C ¥' F (^)- 



Using lemma [5Tl .21 



sm(br])f I,F (x, 77) 

= 7^W + (W) (7 ' C) -no,- b) f- (r (0 ,-6)/ (/lC) )] / ' J? (x,7/) 

= ^[2/ fe e -r ( o,- f ,)(/ + / (/ ' C) )] / ' F (x,r / ) = - r (0 , -^/of'^, 17). 

□ 

This means that the data from two antennas with distance b is sufficient to recover 
sin(br])f I,F . 

5.1.4 Remark 

1. f I,F (x,ri) cannot be recovered with two antennas with distance b for 

2. /^(x,0) = (/ 6 e )^(x,0)V6GR 

5.1.5 Corollary 

Let f G y(H n xR). If f£ is known for all b > 0, then f is completely determined. 
Proof: 

If /q is given, the even part of / is known. The odd part is determined for b > 
by 

7;[.f(x,&) - f(x, -&)] = y sin(6?7) y sin(yr])f(x,y)dydr} 

R 

sin(br])f I ' F (x, rf) ^77 



R 
-7 



V^TT 

R 

-1 



27T 



R 

□ 



5.2 Post-processing formulas 

5.2.1 Definition 

Let / G J^(E n x E). Define 

/i 6 (77) := sin(b7])f hF (x,r]). 
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It turns out that a noise resistant algorithm to recover the odd part of the reflectivity 
function / in addition to the even part is not so easily found due to the many zeros 
of the sine. Some obvious regularizations such as 



(fa) ' 0, 0) for 7] = with a G {0, b} 



sin( y^ otherwise, 

Bin {brfi+e 



(faY' F (x, 0) for r] = with a G {0, b} 



sm(b-q)h b (ri) 
sin 2 (bri)+e cos 2fe (brj) 



otherwise 



or 



(f*Y' F (x,0) for 7] = with a G {0,6} 

sin 2 (b v )+e L% (t^exptv ) otherwise 



with fc,leN and e > yield only high-pass-filtered images. The regularization 

(faY' F (x, 0) for 7] = with a G {0, b}, 



sm(bri)h b (ri) 



sin 2 (b?y)+e cos 2fc (t>77)(exp(rf/y 2i )-l) 



otherwise 



did not provide satisfactory results either. To avoid the problem of a vanishing 
denominator, two satisfactory solutions exist which will be given in the following 
theorem. 

5.2.2 Theorem 

Let f G y(R n x E). 

1. For two data sets, i.e. only one b G R, a regularization is necessary: 

(faY' F (x, 0) for7] = with a G {0, b} 



f' F (x,v) 



sin(br))h b (ri) 

sin 2 (br])+H n (r?)e cos 2 * (677) 



otherwise 



with e > 0, k G N, and the Heaveside function H = X(^-,oo)- 

2. Form > 2 data sets, i.e. a set b\, 6 m ( m -i)/2 £ K, z7 is possible to choose 
the bk in a way that the denominator does not vanish, e.g. b\ = 1, 62 = tt- 



76 



Therefore an analytically exact formula is achievable: 



r (/fj /,J? M) 



m(m— 1) /2 

k=l 

m(m-l)/2 

E sin 2 (6 fe 7?) 

fc=i 



for rj = wzY/z 

G {l,...,m(m- l)/2} 



otherwise. 



Proof: 

1. Let rj £ f b Z, k e N, and e > 0. Then 

sin(6?7)/z fe (ry) 



sin (brj) + HjL(r])ecos 2k (brj) 



sm 2 (bri)f I,F (x, rj) 



f^ F (x, V ). 



sin (677) + HjL(r])ecos 2k (br]) e^o 



2. Let 77 0. Then 



m(m-l)/2 

E sm(b k rj)h bk (r]) 

k=l 

m(m— 1) /2 

E sin 2 (^77) 
fc=i 

m(m— 1)/2 

E sm 2 (b k r])f I ' F (x,7]) 

k=l 

m(m— 1)/2 

E sin 2 (6 fc ?7) 
fc=i 



/^(X,77). 



□ 

5.2.3 Remark 

1. The same as in theorem \5.2.2\ could be shown for the other aforementioned 
regularizations although they do not yield satisfactory results. However this 
property is necessary to be a viable approach. 

2. It is essential not to regularize everywhere. Otherwise only a high-pass 
filtered image is obtained. 

3. Considering the sampling-theorem, the formula for two data sets is exact if 



supp h b C 
than b. 



b ' b 



), i.e. if the bandwidth of f in y -direction is smaller 
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4. An advantage of using more than two data sets is a sizable reduction of 
noise-effects by the weighted summation of the data. 



5.3 Numerical simulations 

In the following, two different kinds of numerical simulations will be shown. First, 
to demonstrate the ability of the proposed algorithms in solving the left-right am- 
biguity, a set of left-right symmetric images will be created that serve as input for 
the algorithms. In this way artifacts arising from the inversion of the spherical 
Radon transform are omitted, and it is possible to isolate the effects of the post- 
processing formulas. Second, to show what results can be expected in reality, 
measurements with multiple antenna positions will be simulated with the spher- 
ical Radon transform that serve as data for reconstructions according to theorem 
12.3.61 Then the acquired left-right symmetric images are fed into the algorithms 
to solve the left-right ambiguity. 

In the following figures antenna distances will be given in multiples of the sam- 
pling length (pixels). The noise is applied as a pixelwise multiplication with the 
corresponding percentage multiplied with a completely uncorrelated (with regard 
to pixels as well as pictures), normally distributed deviate with zero mean and unit 
variance. Additionally a similarly structured absolute error is applied. 
Figure 15.11 shows on the left hand side a simple circle with five smaller circles 
arranged in a cross formation in it. The flight track, i. e. the leftmost antenna, 
runs in a vertical line right in the middle of the figure with the additional antennas 
on the right. In the following this will serve as the phantom for the attempts to 
reconstruct the reflectivity function /. The right hand side of figure I5TT1 shows the 
corresponding cross section. 




Figure 5.1: Reflectivity function / 
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5.3.1 Solving the left-right ambiguity with computed symmet- 
ric images 

In the following examples the even part of / with respect to several flight tracks 
is computed directly. These even images are combined to reconstruct the ground 
reflectivity function / according to the corresponding formulas in theorem T5.2.21 
As an example, figure 15.21 shows the symmetric input data, where the number 
of circle pairs depends on the number of antennas employed and the distance 
between the two circles of the pair is associated with the antenna distance. The 
noise level is reflected in the intensity of the deviations and will be kept at a level 
of 10% in this subsection. 



fS(0,y),fi(0,y),fi(0,y) 




Figure 5.2: The even part of / with respect to three different flight tracks with 
10% noise 



First, the influence of the antenna distance is demonstrated. Figures 15.31 to 15.61 
contrast the achievable results for reconstructions with a small antenna distance to 
reconstructions with a large antenna distance. 

In figure 15.31 the reconstruction for two antennas with a distance of one is de- 
picted with the corresponding image in the Fourier space given in figure 15 .41 The 
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y 



256 



Figure 5.3: Two antennas with distance 1, 10% noise, spatial domain 



x 





Figure 5.4: Two antennas with distance 1, 10% noise, frequency domain. Left: 
Logarithmic scale 



x 
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.20 
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256 

Figure 5.5: Two antennas with distance 9, 10% noise, spatial domain 

7^(0,77) 



x 




Figure 5.6: Two antennas with distance 9, 10% noise, frequency domain. Left: 
Logarithmic scale 
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-128 



x 



phantom is acceptably reconstructed and the mirror image is decently suppressed. 
Unfortunately an error in the low frequencies is annoying that is clearly visible 
from the large values of f I,F around 77 = in the cross section of figure [5~4l and 
the notable slope in the cross section of figure 15731 Figures 1531 and [5761 show the 
reconstruction and its Fourier transform for two antennas with a distance of nine. 
In figure 15751 the phantom seems to be out of focus with shadows to both sides. 
Nevertheless the suppression of the mirror image is satisfactory. As can be seen 
from the steep slope in the cross section in figure [5751 and the large values of f I,F 
around 77 = in the frequency domain in figure I574l in comparison with figure [5761 
small antenna distances lead to problems with the low frequencies. For large an- 
tenna distances there are many frequencies missing on periodical, vertical stripes, 
as becomes evident from figure 15.61 This is an effect of the regularization delin- 
eated in theorem [5.2.2L 1 . An advantage is that there are almost no problems with 
low frequencies. 

This points to the possible benefits of a combination of more than two antennas, 
as it is possible to choose one small and one large distance. The achievable im- 
provements are demonstrated in figures 15771 to 15 .101 with three and four antennas. 

f(0,y) 





256 



-128 



x 



Figure 5.7: Three antennas with positions 0, 2, and 5, 10% noise, spatial domain 



+0 



+- 

^ 2 



±7T 
77 



-0 



Figure 5.8: Three antennas with positions 0, 2, and 5, 10% noise, frequency do- 
main. Left: Logarithmic scale 



In figures [5771 and 15 .81 the reconstruction and its Fourier transform are depicted for 
three antennas with distances two and five. The reconstruction quality improved 
in comparison with the preceding figures, but there is still a slight slope visible in 
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Figure 5.10: Four antennas with positions 0, 1, 4, and 9, 10% noise, frequency 
domain. Left: Logarithmic scale 

the cross section of figure 15.71 However, there are no frequencies missing since 
a regularization is no longer necessary due to data from more than two antennas 
(see theorem r5.2.2l) . In figures I5T91 and 15 .101 the result from four antennas with dis- 
tances one, four, and nine can be seen. The quality of these figures is very good, 
since the phantom is reconstructed very accurately and the mirror image is only 
slightly visible. As can be seen from the preceding images, the reconstruction 
quality improves with more antennas. Since there are problems, if only two an- 
tennas are used, regardless of the distance, it is advisable to employ at least three 
antennas to avoid these problems. The effects of more antennas with regards to 
noise will be further discussed in subsection !5.3.3[ 

5.3.2 Solving the left-right ambiguity with reconstructed sym- 
metric images 

In the following examples measurements for several parallel flight tracks are sim- 
ulated using the spherical Radon transform. Then the corresponding even images 
are formed with the formula in theorem 12.3.61 Finally, the images are combined 
to obtain the ground reflectivity function / using the formulas in theorem r5.2.2[ 
As an example, figure 15.1 II shows the data from the simulated measurements 
where the number of boomerang-shapes depends on the number of antennas em- 
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ployed and the positions of the shapes are related to the antenna distances. The 
noise level will be at a level of 10% throughout this subsection. 
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Figure 5.11: Measured data from three antennas with 10% noise 



In the following, the effect of the antenna distance on the reconstruction quality is 
examined. Figures 15.121 to 15.151 contrast the reconstruction with a small antenna 
distance to a reconstruction with a large antenna distance. 
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Figure 5.12: Two antennas with distance 1, 10% noise, spatial domain 
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Figure 5.13: Two antennas with distance 1, 10% noise, frequency domain. Left: 
Logarithmic scale 
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Figure 5.14: Two antennas with distance 9, 10% noise, spatial domain 
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Figure 5.15: Two antennas with distance 9, 10% noise, frequency domain. Left: 
Logarithmic scale 
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Figures [5.121 and l5 . 1 31 display the reconstruction and its Fourier transform for two 
antennas with a distance of one. The reconstruction of the phantom is satisfactory 
as is the suppression of the mirror image. In comparison with figures 1531 and 15 .41 
which show the corresponding images with the same number and placement of 
the antennas that are reconstructed from computed symmetric images, the error 
in the low frequencies for small antenna distances is negligible, however figure 
15.121 is affected by circular artifacts that arise from the noise in the inversion of 
the spherical Radon transform. The small error in the low frequencies in figure 
15 .121 comes from a suppression of the low frequencies due to limited data in the 
inversion of the spherical Radon transform. Figures 15.141 and 15.151 show the re- 
construction for two antennas with a distance of nine and the corresponding image 
in the frequency domain. The phantom in figure I5TT41 is very blurry - an effect of 
the many missing frequencies that can be observed in figure 15.151 Surprisingly 
however, the amplitude of the phantom is more exactly reconstructed in compar- 
ison to figure 15. 121 This points again to a more exact reconstruction of the low 
frequencies for large antenna distances, although the effect is not as striking as in 
subsection 15.3.11 This is also reflected in the slightly smaller maxima in figure 
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Figure 5.16: Three antennas with positions 0, 2, and 5, 10% noise, spatial domain 
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Figure 5.17: Three antennas with positions 0, 2, and 5, 10% noise, frequency 
domain. Left: Logarithmic scale 



Again, this emphasizes that it could be an interesting idea to combine more than 
two antennas as it could allow to gain the advantages of a small and a large antenna 
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distance. The successes are demonstrated in figures I5TT61 to 15 . 1 91 
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Figure5.18: Four antennas with positions 0, l,4,and9, 10% noise, spatial domain 
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Figure 5.19: Four antennas with positions 0, 1, 4, and 9, 10% noise, frequency 
domain. Left: Logarithmic scale 



In figures l5TT6l and l5TT7l the results from three antennas with distances two and five 
can be seen. The reconstruction is superior to the preceding two, but the amplitude 
of the phantom is still not completely correct, as can be seen in the cross section 
of figure 15.161 Again, with three antennas there are no frequencies missing, as 
seen in figure I5TT71 since no regularization is necessary. Figures [5.181 and f5. 191 
show the reconstruction and its Fourier transform for four antennas with distances 
one, four, and nine. The reconstruction improved over the preceding one, but still 
the amplitude of the phantom is not correct. This is an effect of the inversion 
of the spherical Radon transform with limited data. Nevertheless the suppression 
of the mirror image is good. As could be presumed, the low frequencies are 
reconstructed quite nicely and there are no frequencies missing, as can be seen in 
all four figures. The reconstruction quality improves further with more antennas, 
as can be seen in comparison of figures l5T2ll5.14l[5.16l and !5.18l This is not only 
an effect of averaging out the noise, but is amplified by a weighted summation of 
the data, where the data is more heavily weighted the more noise resistant it is. 
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5.3.3 Analysis of stability and the effects of noise 



Now an analysis of noise effects and stability is performed. 
Figures 15.201 to 15.251 show reconstructions from two antennas with varying 
amounts of noise. As in subsection 15 .3 . 1 1 figures 15.201 to 15.221 demonstrate only 
the characteristics of the solution of the left-right ambiguity, whereas figures [5T231 
to [525] contain also the artifacts arising due to the inversion of the spherical Radon 
transform with limited data as in subsection l5.3.2l 
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Figure 5.20: Two antennas with positions and 1, no spherical Radon transform, 
10% noise 
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Figure 5.21: Two antennas with positions and 1, no spherical Radon transform, 
20% noise 
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Figure 5.22: Two antennas with positions and 1, no spherical Radon transform, 
30% noise 



Figures 15.201 to 15.221 show the reconstructions from computed even images with 
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respect to two antennas with a distance of one. In figure 15301 the noise level is at 
10%, in figure EH at 20%, and in figure [5221 at 30%. In figure E^Ql the phantom 
is clearly visible despite the known error in the low frequencies. The phantom in 
figure 15.211 is obscured more strongly due to this problem, and in figure 15.221 it 
is very difficult to discern the phantom from the underlying error. Nevertheless 
it can be seen in the cross sections that the algorithm is stable since the cross 
sections look quite similar disregarding the scaling, and there are no signs of an 
overamplification of noise. 
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Figure 5.23: Two antennas with positions and 1, with spherical Radon transform, 
10% noise 
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Figure 5.24: Two antennas with positions and 1, with spherical Radon transform, 
20% noise 
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Figure 5.25: Two antennas with positions and 1, with spherical Radon transform, 
30% noise 



In figures 15.231 to 15.251 the reconstructions using the inversion of the spherical 
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Radon transform for two antennas with a distance of one are depicted. In figure 
[?T23l the noise level is at 10%, in figure[534]at 20%, and in figure[535]at 30%. The 
reconstruction quality in figure 15.231 is again acceptable, but declines for higher 
noise levels as in figure 15341 Finally, for a noise level of 30% the reconstruction 
quality is not very good, but the phantom is still discernible in contrast to figure 
15.221 All three images do not display the right amplitude of the phantom, but 
as previously analyzed, this is an effect of the inversion with limited data. With 
increasing noise levels the reconstruction degrades, but no uncontrolled noise am- 
plification can be observed. 

Figures [5 .201 to [5T251 show that the reconstruction quality gets worse with a higher 
amount of noise, since the amplitude of the error gets larger. Nevertheless there 
are no signs of an overamplification of noise. The loss of reconstruction quality 
is within the expected order of magnitude, so this points to the stability of the 
algorithm. 

Figures l5.26l to l573"Tl compare the results for calculations with three antennas with 
different noise severity. Figures 15.261 to [5T28l are directly computed from symmet- 
ric images, whereas figures l572~9~l to l5.3 fl are gained from inversions of the spherical 
Radon transform. 

Figures 15 .261 to [5T281 show the reconstructions using computed symmetric images 
with regard to three antennas with distances one and three. The noise level in 
figure [5T261 is at 10%, in figure [3371 at 20%, and in figure [5T281 at 30%. The 
reconstruction of the phantom in figure [5T261 is satisfactory since the amplitude is 
close to the correct values. The suppression of the mirror image is also satisfying. 
However there is still a problem with low frequencies as can be seen in the slight 
slope in the cross section in figure 15361 The reconstruction quality in figure [5371 
is worse. The phantom has a higher peak and the mirror image is more strongly 
visible. Additionally, the slope in the cross section of figure [5371 is steeper which 
points to a more severe problem with low frequencies. This tendency continues in 
figure 15381 The amplitude of the phantom is again farther off, albeit not by much, 
and the mirror image is more pronounced. Also the slope in the cross section of 
figure 15 .281 steepened a bit. 

In figures 15.291 to 15.311 the results for calculations with three antennas with dis- 
tances of one and three using the inversion of the spherical Radon transform are 
displayed. In figure [5391 the noise level is at 10%, in figure [57301 at 20%, and in 
figure I573T1 at 30%. The reconstruction in figure [5391 is quite good. The ampli- 
tude of the phantom is not completely correct, but as already mentioned this is an 
effect of the inversion of the spherical Radon transform due to limited data. The 
suppression of the mirror image is also acceptable. In figure [57301 a clearly inferior 
reconstruction quality due to notable errors in the amplitude of the phantom and 
less suppression of the mirror image can be observed. Figure I573T1 is similar to 
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Figure 5.26: Three antennas with positions 0, 1, and 3, no spherical Radon trans- 
form, 10% noise 
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Figure 5.27: Three antennas with positions 0, 1, and 3, no spherical Radon trans- 
form, 20% noise 




Figure 5.28: Three antennas with positions 0, 1, and 3, no spherical Radon trans- 
form, 30% noise 
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Figure 5.29: Three antennas with positions 0, 1, and 3, with spherical Radon 
transform, 10% noise 




Figure 5.31: Three antennas with positions 0, 1, and 3, with spherical Radon 
transform, 30% noise 
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figure 15.301 but overall the quality is a bit worse due to higher variances in the 
phantom's amplitude and the mirror image. 

Figures 15.261 to l53T1 show comparable results to figures [5T20l to 15.251 but the qual- 
ity is clearly better. The images with 10% noise are satisfactorily reconstructed. 

In order to examine, whether the reconstruction quality - especially with respect 
to noise - can be enhanced with even more than three antennas, in the following 
the results of calculations with four antennas will be discussed. 
Figures l5.32l to l5371 dispray reconstructions from four antennas with several noise 
amplitudes. Figures [5.321 to [5341 show only the effect of the solution of the left- 
right ambiguity whereas figures 15.351 to 15.371 also include the artifacts from the 
inversion of the spherical Radon transform. 

Figures [5.321 to [5341 display reconstructions from computed even images for four 
antennas with distances one, three, and eight. The noise level in figure I532l is at a 
level of 10%, in figure [5331 at 20%, and in figure [5341 at 30%. The reconstruction 
in figure [5321 is very good. The amplitude of the phantom is recovered very well 
and the mirror image is quite small. Only a very slight slope can be detected in the 
cross section. The quality in figure [5331 declines, but nevertheless the amplitude 
is almost on target and the suppression of the mirror image is still acceptable. 
Noticable however is the increased slope in comparison to figure 15.321 Figure 
15.341 shows an additional decline in quality, but not as severe. The amplitude 
deviates farther, and the slope is steeper. Most annoying however is the more 
strongly pronounced mirror image. 

Figures [535] to [537] delineate the results from the inverted spherical Radon trans- 
forms using four antennas with distances one, three, and eight. In figure [5351 the 
noise level is at 10%, in figure [5361 at 20%, and in figure [5371 at 30%. Figure 
15.351 looks quite good. The suppression of the mirror image is adequate and the 
amplitude of the reconstructed phantom is almost only plagued by the error due 
to the inversion of the spherical Radon transform with limited data. In figure [5361 
the suppression is not as good and the amplitude of the phantom deviates farther. 
Finally figure [5371 displays a slightly larger deviation in the object's amplitude 
and a more noticable mirror image. 

Figures 15.321 to 15.371 show a sizeable improvement with regard to figures 15.261 to 
I5.311 Now, even with 20% noise the reconstructions seem quite useful and the 
images with 10% noise look rather good. 

But this effect can be enhanced even more. 

Figures [5.381 to I5~43l are reconstructions from five antennas with a maximum dis- 
tance of 19 and varying amounts of noise. With regard to currently achievable 
resolutions and the wingspan of the used airplanes, this distance corresponds to 
an appropriate maximum antenna distance. Figures 15.381 to 15.401 are computed 
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Figure 5.32: Four antennas with positions 0, 1, 3, and 8, no spherical Radon 
transform, 10% noise 
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Figure 5.33: Four antennas with positions 0, 1, 3, and 8, no spherical Radon 
transform, 20% noise 
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Figure 5.34: Four antennas with positions 0, 1, 3, and 8, no spherical Radon 
transform, 30% noise 
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Figure 5.35: Four antennas with positions 0, 1, 3, and 8, with spherical Radon 
transform, 10% noise 




Figure 5.37: Four antennas with positions 0, 1, 3, and 8, with spherical Radon 
transform, 30% noise 
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from symmetric images. Figures 15 .4 II to [5~43l are the results of a complete recon- 
struction process, including the inversion of the spherical Radon transform. 
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Figure 5.38: Five antennas with positions 0, 1,3, 8, and 19, no spherical Radon 
transform, 10% noise 
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Figure 5.39: Five antennas with positions 0, 1,3, 8, and 19, no spherical Radon 
transform, 20% noise 
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Figure 5.40: Five antennas with positions 0, 1,3, 8, and 19, no spherical Radon 
transform, 30% noise 



Figures 15.381 to 15.401 show reconstructions from computed even images with re- 
spect to five antennas with distances one, three, eight, and 19. The noise level in 
figure 15381 is at 10%, in figure 15391 at 20%, and in figure [5391 at 30%. Figure 
15381 exhibits an almost perfect reconstruction. The object's amplitude is recov- 
ered very accurately, and no slope is visible. The most notable feature however is 
that the mirror image is almost undiscernible from the noise floor. Also in figure 
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15.391 a very good reconstruction can be seen. The amplitude of the phantom is 
quite good, and again there is no slope. The mirror image is visible, but does not 
stand out excessively. The reconstruction quality visible in figure [5401 is still re- 
markable. The amplitude of the object is recovered closely and there is no slope. 
The suppression of the mirror image is also very good. 
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Figure 5.41: Five antennas with positions 0, 1, 3, 8, and 19, with spherical Radon 
transform, 10% noise 
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Figure 5.42: Five antennas with positions 0, 1, 3, 8, and 19, with spherical Radon 
transform, 20% noise 
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Figure 5.43: Five antennas with positions 0, 1,3, 8, and 19, with spherical Radon 
transform, 30% noise 



Figures 15.411 to 15.431 delineate reconstructions using inversions of the spherical 
Radon transform from five antennas with distances one, three, eight, and 19. In 
figure |5~4T1 the noise level is at 10%, in figure [5421 at 20%, and in figure [5431 at 
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30%. The reconstruction depicted in figure [5~4T1 is very good. The mirror image 
almost drowns in the noise induced by the inversion of the spherical Radon trans- 
form with limited data. Unfortunately there is still a gradient in the amplitude of 
the phantom. However, as can be seen in comparison with figure 15381 this is ob- 
viously an effect of the limited data inversion, too. Figure 15421 is still plagued by 
the problems of the limited data inversion, but apart from that the reconstruction is 
quite good. The suppression of the mirror image is still quite good. However the 
phantom's amplitude varies too much. Figure 1543] displays a stronger deviance in 
the object's amplitude. The mirror image is also more pronounced. Therefore the 
reconstruction quality is a bit worse. 

Figures 15 .3 81 and 15 .4 II show only very scant signs of noise. But also figures [5391 
and 15.421 are acceptable and figures 15.401 and 15.431 show a large improvement in 
comparison with reconstructions that use less antennas. Figures l5.38l to l5431 show 
a noticable improvement in comparison with figures [5321 to 15.371 The most im- 
portant improvement is the absence of a slope in figures [5381 to 15.431 but also the 
better suppression of the mirror image is obvious in all figures. 
The results of subsection [533] show that the algorithms in theorem [533] are sta- 
ble. Additionally, they show that it is advantageous to employ more than two 
antennas and that better results can be expected for an increasing number of an- 
tennas. So if an airplane would be used to its maximal capability, then even with 
a large interference by noise very good images should be possible and the results 
should be very stable. This culminates in the insight that even for very adverse 
noise levels it is possible to receive good images if the whole wingspan of an 
airplane can be used. 
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Chapter 6 

Conclusions and outlook 



In this thesis several topics connected with SAR are addressed. Most of them will 
help to increase image quality, but one raises questions that deserve further study. 
In chapter [2] an important error in a common approach to invert the spherical 
Radon transform was analyzed and a solution was proposed. This allowed for a 
promising new approach to alleviate the problem of limited data. The results of 
numerical simulations are very promising. To improve the results even further, it 
should be studied whether and what kind of decay in the approximatively contin- 
ued data enhances the reconstruction quality. 

Chapter[3]addressed the problem of limited data more thoroughly. The findings of 
this chapter contribute to an understanding of the artifacts that frequently appear in 
the numerical inversion of the spherical Radon transform. Unfortunately, it turned 
out that the artifacts caused by limited data are not restricted to high frequencies, 
as is the case in computerized tomography. Since this analysis should also hold 
for other models of SAR, e.g. models using the wave equation, the implications 
of this finding should be examined with great care. As can be seen in the numer- 
ical simulations of chapter [2j the applied regularization has a strong influence on 
the artifacts that appear in the reconstruction. Therefore the physical meaning of 
regularizations should be studied. Hopefully, this results in a reconstruction for- 
mula that minimizes artifacts and that can be relied upon, because it matches the 
physical conditions. 

Chapter @] presented another new way to invert the spherical Radon transform 
that uses the ideas of chapter [3] to reduce the complexity of the problem with a 
projection onto orthogonal functions. With the insight gained from chapter [3] it 
might be possible to select some of these orthogonal functions to obtain good 
reconstructions. If only a handful of these function would be accounted for, this 
could also yield a fast way to invert the spherical Radon transform. 
Chapter [5] dealt with the problem of the left-right ambiguity that causes objects 
that lie on one side of the flight track to appear in the reconstruction of both sides. 
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With the postprocessing formula derived in this chapter it is possible to use recon- 
structions from at least two antennas to solve this problem. The numerical sim- 
ulations are very promising. It can be concluded that the reconstruction quality 
improves with a higher number of antennas. It is however important to note that 
this is not only the usual effect of averaging over more data to reduce noise. This 
effect is amplified by a weighted summation of the data that emphasizes data less 
affected by noise. Therefore if an ample number of antennas would be mounted 
on an airplane, it should be possible to reconstruct images that show only very 
scant errors despite heavy noise. 

Because of the jitter of the airplane a modification of the formulas in this chap- 
ter should be studied. It should be possible to obtain a similar formula that uses 
the data instead of the reconstructed images as input. However it would be nec- 
essary to mirror the data so it becomes an even function that has a meaning for 
negative radii. Then the left-right ambiguity in the data could be treated in a sim- 
ilar way as in chapter [51 and for each side of the flight track a separate inversion 
could be performed. With this modification the fluctuation in airplane position 
and heading would probably no longer pose a problem because it is very small 
for one send/receive cycle. An additional advantage would be the fact, that the in- 
version of the spherical Radon transform would have to be performed only twice 
regardless of the number of antennas. This should speed up the reconstruction 
considerably. 
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